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EXECUTIVE SUMMARY 


The results of studies of wind shear hazards to aircraft operation 
carried out under NASA Marshall Space Flight Center contract for the 
period 1979 through 1981 are summarized in this report. The results of 
the study are integrated with other reported information in the litera- 
ture and with cooperative programs carried out with NASA Ames Research 
Center and United Airlines Flight Training Center. 

The report first reviews existing wind shear profiles currently 
used in computer and manned flight simulator studies. The governing 
equations of motion for an aircraft are then derived incorporating the 
variable wind effects. Quantitative discussions of the effects of wind 
shear on aircraft performance are presented. These are followed by a 
review of mathematical solutions to both the linear and the nonlinear 
form of the governing equations. Solutions with and without control 
laws are presented. 

The application of detailed analysis to developing a warning and 
detection system based on a Doppler radar measuring wind speeds along 
the flight path is given. These real-time wind speed profiles are fed 
into a microcomputer, and utilizing the governing equations of aircraft 
motion, a flight path deterioration parameter representing a measure of 
the severity of the wind shear is predicted. A number of flight path 
deterioration parameters are defined and evaluated. Comparison of 
computer-predicted flight paths with those measured in a manned flight 
simulator for flight through hypothetical sinusoidal wind shears and 1 - 
cosine downdrafts is made. The fidelity of the computer program calcu- 
lations with the measured manned flight simulator aircraft response is 
described. Also a correlation of the magnitude of the flight path 
deterioration parameters with aircraft controllabil ity along the flight 
path for varying magnitudes of sinusoidal wind speed amplitudes and 
frequency oscillations is given. 
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The report ends with a review of some proposed airborne and ground- 
based wind shear hazard warning and detection systems. The advantages 
and disadvantages of both types of systems are discussed. 

The conclusions of the review are that existing wind shear models 
used in computer and manned flight simulator studies are not realistic. 
All existing mathematical models of wind shear are spatially two- 
dimensional and based on highly smoothed and limited data; none include 
time dependence. Moreover, the small-scale microburst-type wind shear 
is not contained in any of the models. Complete data sets from which 
very good wind shear models can be developed are now available through 
the NASA Gust Gradient and NCAR JAWS field programs, but these need to 
be analyzed. 

Order of magnitude analysis of the equations of motion for an 
aircraft illustrates that low values of horizontal wind shear are much 
more hazardous than larger values of vertical wind shear. The FAA 
AC-20-57A Advisory Circular , relative to the certification of automatic 
control systems, calls for 8 kts/100 ft but does not specify that the 
value be measured along the flight path. The value implies 8 kts/100 ft 
of altitude. It is believed that realistic three-dimensional time- 
dependent wind shear models should be used for certification. 

Argument exists as to the correct flight procedure to employ when 
caught in severe wind shear. The main controversy is relative to the 
optimum speed to fly during an encounter with a head wind shearing to a 
tail wind. Controversy as to whether to fly at stick-shaker speed or 
minimum drag speed exists. The Aline Pilots Association (ALPA) Air- 
worthiness and Performance Committee recommends flying at minimum drag 
speed and thus maintaining some excess kinetic energy to flair the 
aircraft at the last instant if impact cannot be avoided. 

Initial calculations of flight through wind shear showed conflict- 
ing results depending on whether the wind speed profile varied linearly 
in the vertical or logarithmically. This disagreement can be traced to 
the initial or trimmed condition used in the analysis. 
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Computer simulation of aircraft flying through several mathematical 
thunderstorm models developed from gust front data measured with a 500-m 
(1500 ft) tower at NOAA/NSSL, Norman, Oklahoma, clearly illustrates that 
the amplitude of the phugoid oscillation of the aircraft is highly 
amplified. Small perturbation stability analysis clearly supports this 
observation. Because the wind shear in thunderstorms creates a force 
function having essentially the same frequency of the aircraft phugoid, 
it is believed that the phugoid mode normally considered benign can 
become hazardous when flying through a thunderstorm. Careful evaluation 
of flight training simulators to assure valid reproduction of the air- 
craft phugoid characteristics should be made when using simulators to 
train flight crews or evaluate airborne systems. 

Flight path deterioration parameters computed from wind speeds 
measured with the Doppler radar looking along the flight path show good 
promise as an effective index of hazard level for use in wind shear 
warning and detection systems. Comparison of flight path deterioration 
parameters evaluated through computer simulation with those measured 
through manned flight simulators (i.e., with man in the loop) show 
generally consistent results. Additional work is required, however, to 
establish a meaningful magnitude of the parameter and a scale of wind 
shear severity. 

Airborne wind shear warning and detection systems have been evalua- 
ted and have proven effective in manned flight simulator studies. The 
airborne aids, however, have been tested primarily for approach flight 
conditions using the standard wind shear models which are not believed 
to represent realistic nor the most severe conditions. The basic 
principle of the airborne aid is to maintain ground speed thus storing 
energy for conditions when the head wind shears to a tail wind. This 
system, of course, has limited use during takeoff at essentially maximum 
thrust. Additionally, the airborne system has the disadvantage that one 
must be in the wind shear before it provides any warning. Finally, the 
need for a very accurate ground speed measurement, not normally avail- 
able on board the aircraft, is required for these systems. 
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The low-level wind shear alert system, LLWSAS, which is a ground- 
based warning and detection device and which has been installed at 
58 major airports as of October 1982, must be considered only an interim 
solution. Current studies have clearly indicated that the scales of 
extreme wind shear are sufficiently small such that they can go undetec- 
ted by most LLWSAS arrangements. Additionally, these are a surface 
measurement and do not provide warning when wind shear occurs along the 
flight path but outside the airport perimeter. Finally, arguments have 
been made that they give too many false alarms resulting in their 
warning being ignored in many cases. 

The opmum warning system appears to be a dual Doppler radar which 
has been demonstrated without reservation to be capable of monitoring 
all necessary scales of wind shear. The cost of installing Doppler 
radar at every major airport may be prohibitive. 

It has not yet been resolved as to whether monitoring the component 
of wind along the flight path, which all current ground-based and 
airborne detection systems do, is adequate. The vertical component of 
the wind may be a very significant parameter which must also be moni- 
tored. This, of course, can be measured using two Doppler radars; 
however, the cost of installation is compound. Further study as to 
whether the effect of the downdraft on airplane performance can be 
ascertained by monitoring only the longitudinal wind speed component 
and as to meaningful magnitudes of the downdraft velocity close to the 
ground is needed. 



1.0 INTRODUCTION 


With the advent of the digital flight data recorder, low-level wind 
shear has been recognized as a severe flight hazard [1]. Investigations 
of at least 25 commercial airline accidents and at least 5 U.S. Air 
Force (USAF) mishaps [2,3] have clearly proven that wind shear, result- 
ing in a sudden change in either the speed or direction of the wind, can 
produce dynamic effects on aircraft which cause them to deviate signi- 
ficantly from the pilot's intended flight path producing impact with the 
ground or frightening near-misses. Both the International Civil Aviation 
Organization (ICAO) and the Federal Aviation Administration (FAA) now 
recognize wind shear as a potential hazard to the safety of aircraft 
operations, especially in the critical landing and takeoff phases of 
flight. Prior to this recognition, the role wind shear played in 
aircraft accidents may often have been attributed to pilot error. 

It is not surprising that the temporary loss of control or struc- 
tural failure due to unusual and extreme wind variations has gone 
undetected for many years. Practically all textbooks (see as examples 
References 4 through 8) and education programs on aircraft flight 
dynamics consider only constant or zero winds both in the development of 
the governing equations and in the analyses of aircraft motion in the 
atmosphere. It should be noted that although numerous studies relative 
to the influence of individual gusts or random turbulence on flight 
performance of aircraft (see for example References 9 through 13) have 
been conducted, these are generally associated with the high-frequency 
atmospheric fluctuations. Thus, only aircraft performance relative to 
changes in wind on time and spatial scales, which are small in compari- 
son with the scales associated with severe wind shear (see for example 
References 14 and 15), have been studied. Moreover, only recently is 
wind shear of this scale being measured in the detail necessary to 
analyze its effect on the motion of aircraft [16,17]. 
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still, however, insufficient meteorological data are available to 
construct three-component, three-dimensional spatial, and time-dependent 
models of wind shear for aircraft design, operational procedures develop- 
ment, and simulation studies. Models of wind fields associated with 
thunderstorms and other sources of extreme wind shear are urgently 
needed to develop and verify existing detection and warning systems, to 
upgrade manned flight simulators for training purposes, and to establish 
structural and control design criteria. 

The purpose of this report is to document and compile information 
on aircraft procedures and safety during operation in a wind shear 
environment. Much of the information was developed and assembled under 
NASA-supported or jointly-supported programs. 

The report first describes existing wind shear models and indicates 
where additional data are needed. Next it summarizes some of the 
effects on aircraft performance due to spatial and temporal variation in 
the wind. The dynamic equations of motion are developed and additional 
terms, which occur due to the variable wind effects, are described. 

Some simple calculations are made to illustrate the influence of these 
additional terms on typical approach and takeoff through wind shear. A 
review is then given of previous studies of the effect of wind shear on 
aircraft performance. In these studies, a number of restrictive assump- 
tions, such as linearity of the wind shear profile, variation only in 
horizontal winds, or three-degrees-of-freedom motion are made. Results 
of analyses of more extreme wind shear conditions, such as have been 
associated with aircraft accidents, are then reviewed. Finally, recent 
studies relative to the development of detection and warning systems are 
described. 
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2.0 WIND SHEAR MODELS 


2.1 Needs for Improved Wind Shear Data 

The need for additional wind shear data is manifest in manned 
flight simulator studies, structural and control design analyses, and 
detection and warning systems development. 

The FAA [18] proposes to permit expanding training, checking, and 
certification of flight crew members in advanced flight training simu- 
lators. Under the advanced simulation plan, the simulators will have 
the capability to be programmed to represent a full range of aircraft 
flight conditions, as well as specific aircraft accidents in abnormal 
environmental conditions. In this way, flight crews can experience a 
far-ranging set of flight environments and malfunctions, which will 
assist the crew in making proper judgments when abnormal situations 
occur in flight. 

Phase II of the FAA proposed simulator upgrade program includes the 
requirement for representative crosswind and three-dimensional wind 
shear dynamics based on aircraft-related data. In another FAA report 
[19], seven candidate standard wind shear profiles for systems quali- 
fications are reported. These model s--al though fast becoming standards — 
are not truly three-dimensional wind shear. They were constructed from 
data measured with instrumented towers, from reconstruction of winds 
from accident flight data records, and from meteorological math models. 

2.2 Current Wind Shear Models 

The proposed seven candidate standard wind models were selected 
from 21 models investigated with computer and manned flight simulator 
studies [19]. They consist of one mathematical model, three tower 
measurements, and three accident reconstructions. 
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The mathematical model is considered a low-severity wind shear 
condition and represents neutral atmospheric conditions. The three 
tower measurements, one from Cedar Hill tower data [20] and two from the 
500-m tower at the National Severe Storms Laboratory (NSSL), Norman, 
Oklahoma [17], are considered to be of low to moderate severity. The 
tower data from Cedar Hill are considered to represent a nighttime 
stable boundary layer, whereas the data from NSSL represents thunder- 
storm conditions. The three accident reconstructions are the Logan 
International Airport, Boston, 1973, Iberian DC-10 Airline accident; the 
Kennedy International Airport, New York, 1975, Eastern B-727 Airline 
accident; and the Philadelphia International Airport, Philadelphia, 

1976, Allegany DC-9 Airline accident. The latter two accidents are 
considered to represent high-severity thunderstorm models, whereas the 
former is considered to represent winds associated with a warm front of 
moderate severity. 

Although Reference 19 concludes that a collection of realistic 
three-dimensional wind models of three levels of severity have been 
established, variation of the wind field in a lateral direction from the 
flight path and with time is not included. Thus, during a simulation 
with these wind models an aircraft moving sideways to the wind field 
would experience uniform winds in that direction, which is a highly 
unlikely situation. 

Recently, there has been growing evidence that a small-scale but 
severe low-level thunderstorm wind, now referred to as a "microburst," 
occurs with surprising frequency, and cannot only adversely affect 
airplanes but can produce major damage to property on the ground [21]. 
The precise nature of these small-scale events is not clear, but air- 
craft accident investigations and surface damage surveys indicate their 
horizontal extent is typically less than 5 km (3.2 mi) in length and 1 
to 5 min in duration. Unfortunately, most previous thunderstorm inves- 
tigations have not concentrated on such a small scale but rather on the 
larger scale (5 to 25 km), which is more closely related to gust fronts, 
tornado cyclones, and overall storm structure. Because the proposed 
standard wind shear profiles are highly idealized and/or heavily 
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smoothed, it is believed that they do not include the detailed kinematic 
structure of these events. 

To understand the nature of thunderstorm wind shear (probably the 
most hazardous wind shear condition) and the limitations of the current 
models, a description of the thunderstorm is necessary. Figure 2.1 is 
a simplified cross section of a thunderstorm. General airflow and 
precipitation features are indicated. Of particular interest are the 
occurrences of downdrafts and outflow regions, which account for rapidly 
varying winds, or wind shear, in the low levels. Substantial insight 
into the larger scale nature of extended thunderstorm outflow has been 
given by Goff [17], Frost and Camp [23], and Goff [24] in several 
examinations of gust fronts. An expanded view of the outflow or gust 
front region of Figure 2.1 is given in Figure 2.2. Goff et al . [22] 
based his gust front description on measurements of winds during the 
passage of thunderstorms by a 500-m (1500 ft) tower. 

Frost et al . [14] utilized the data from Goff [17] to construct 
tabulated thunderstorm wind fields for use with computer "lookup" rou- 
tines. The data set from Goff [17] consists of longitudinal, W^, 
lateral, W^, and vertical, W^, wind speed components in a vertical 
plane. Data from 20 thunderstorms were measured during the months of 
May through June over the period of 1971 through 1973 with the WKY- 
TV/NSSL 500-m meteorological tower, Norman, Oklahoma. Time histories of 
the wind speeds were converted to horizontal spatial distributions using 
Taylor's hypothesis (i.e., x = Wt) . Ten-second averaged values of wind 
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speed components are provided in the form of isotach maps for W^, W^, 
and W^, respectively. These data were interpolated onto a 41 x 11 point 
grid system, as illustrated for the horizontal wind component in Figure 
2.3, and stored as discrete values on magnetic tape. It should be 
noted, however, as shown by the insert on Figure 2.3, that the data 
represent only a vertical plane through that part of the storm which 
passes the tower. Tabulated values for all 20 thunderstorm wind speeds 
are given in Frost, et al. [16]. The thunderstorm tower data discussed 
in Foy [19] are similarly tabulated. 

Many thunderstorms may not contain well-defined gust fronts (regions 
of outflow extending over many kilometers) as defined by Goff. However, 
essentially all thunderstorms contain downdraft air, which usually 
impacts and spreads out over the surface. Figure 2.4 shows a schematic 
view of a high-plains or desert thunderstorm where dry, cold air tends 
to produce significant downdrafts. Although extensive evidence is 
lacking, in some cases, the downdraft and immediate outflow associated 
with it at the surface can be quite intense and can occur on a rather 
small scale. 

Fujita and Wakimoto [25] and Fujita and Caracena [26] have per- 
formed several analyses on a phenomenon they have termed microburst, 
to indicate a coupled small but intense downdraft and outflow, which 
occurs in thunderstorms which, in most cases, may be of very low 
rainfall rate/radar reflectivity. One such analysis depicts a micro- 
burst occurring along the flight path of Eastern Airlines Flight 66, 
which crashed short of the runway at New York's JFK Airport in 1975. 

This conceptual analysis, using sparse data from the on-board flight 
data recorder, is presented in Figure 2.5. Keenan [27] developed one of 
the proposed standard wind shear models reported in Reference 19 by 
laying a grid system on Figure 2.5. By sketching in flow lines using 
the numbers determined from the flight data recorder and employing 
conservation of mass, a spatial model of the wind field from the Eastern 
66 accident was reconstructed. In this case, both the lateral variation 
of the wind field as well as the lateral velocity component itself are 
unknown. Thus, during a simulation, the aircraft experiences no real- 
istic lateral wind component, and if displaced laterally, it "sees" no 
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Horizontal Grid Points 

Figure 2.3 Grid system superimposed on a typical thunderstorm Wind speed contour map. 
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Figure 2.5 The path of Eastern 66 on June 24, 1975, in the vertical 

plane including the glide slope of runway 22-L at JFK [26]. 



variation in the wind. Insufficient data are available to fully deter- 
mine how strong the wind shear in the lateral direction can be and 
therefore how significantly it can influence an aircraft during approach 
and takeoff. 

Fujita [28] attempted a detailed examination of microbursts in 
Project NIMROD (Northern Illinois Meteorological Research on Down- 
bursts). On May 29, 1978, an interesting observation of an intense 
microburst, which occurred near a National Center for Atmospheric 
Research (NCAR) Doppler radar installation at Yorkville, Illinois, was 
made. Figure 2.6 (taken from Fujita and Wakimoto [25]) shows analyzed 
Doppler velocity fields for this event. The maximum horizontal wind 
measured was 31 m s~^ (60 kts), at an altitude less than 200 m and 
probably as low as 20 to 60 m (66 to 196 ft) above the radar. Such an 
intense microburst occurring so low to the surface would be extremely 
hazardous should an aircraft encounter it during takeoff or approach. 
Other downburst events are reported in Fujita and Wakimoto [25]. 

Two 3° glide slopes are drawn on Figure 2.6. The approximate wind 
speeds along paths #1 and #2 are compared with values reconstructed from 
the Eastern 66 accident (from Foy [19]) and from the tower data. 
Thunderstorm #9 (Case H, Goff [17]) in Figure 2.7. In general, the 
longitudinal wind speed profile is similar but the magnitude is larger 
for the microburst. The very high vertical wind speed proposed by Foy 
[19] is not apparent in either the microburst or the gust front. In 
fact, no other apparent sources of wind shear data show such a strong 
downdraft. This high value is undoubtedly due to assuming only two- 
dimensional conservation of mass when reconstructing the wind field from 
Figure 2.5. 

Frost et al . [14] investigated the magnitude of vertical down- 
drafts. The smoothed data from the 500-m (1500 ft) tower gave downdraft 
values no greater than 3 m/s. Actual values of downdraft as high as 
15.5 m/s (not those reconstructed by Foy [19]) were reported in Fujita 
and Caracena [26]. These values, however, are undoubtedly averaged over 
much shorter periods of time than 10 seconds for which the data pre- 
sented by Goff are averaged. Reference 26 gives no information on the 
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U-COMPONENT WINDSPEED and GROUND-RELATIVE FLOW 



0 I 2 3 km 


Figure 2.6 A vertical cross section through the May 29, 1978, microburst 
showing isotachs of horizontal wind speeds. The height of the 
maximum wind is estimated to be 50 m or lower. Arrows are 
ground-relative velocities in the plane which is stretched 
vertically [25]. (The flight path lines #1 and #2 are super- 
imposed by the present author.) 
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averaging time utilized in arriving at the quoted value of 15.5 m/s. It 
is apparent, however, that the thunderstorm wind field developed by 
Keenan [27] and reported in Foy [19] contain much more extreme down- 
drafts than those measured by Goff [17]. Thus, there are conflicting 
data and opinions as to the maximum magnitude of the downdraft that can 
occur in a thunderstorm. Although it is expected that the 10-sec 
averaged data of Goff will have lower values than the peak downdraft 
wind speed reported by Fujita and Caracena, the discrepancies in the 
values cannot be completely attributed to averaging time. 

Alexander [29] gives a statistical summary of vertical wind speed 
data recorded with NASA's 150-m (500 ft) ground wind tower facility, 
Kennedy Space Center, Florida. One year of continuous around-the-clock 
vertical wind speed measurements were processed to determine the inten- 
sity, frequency, time of occurrence, etc. of the daily maximum vertical 
gust. Both updrafts and downdrafts were studied. These values repre- 
sent 0.1-sec averages, and the maximum vertical downdraft recorded was 
9.3 m/s (18.1 kts), although data recorded during Hurricane Agnes did 
contain an extreme downdraft in excess of 11.9 m/s (23.1 kts). 

Sinclair [30] indicates downdrafts at an altitude of 100 m (300 ft) 
for an Oklahoma thunderstorm may be considerable in excess of 15.5 m/s 
(30 kts). Sinclair has measured downdrafts as high as 28 m/s based on a 
1/25- sec averaging period. Finally, the numerical models of Williamson 
et al. [31] do not predict wind speed downdrafts greater than 10 m/s 
(19.4 kts). Thus, it is evident that research is needed to resolve the 
magnitude of the maximum downdrafts that can occur in a thunderstorm and 
the heights at which they occur. 

2.3 Scales of Wind Shear 

A critical aspect of shear is the length and time scales over which 
the wind is measured. In the atmosphere, with particular reference to 
aircraft problems, four scale regimes can be defined [32]: 

1. The small or turbulent scale which may extend to scales 
as large as a few hundred meters. 
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2. The cloud scale which may range from 0.01 to 10 km 
(0.006 to 6.2 mi ) , 

3. The overall thunderstorm which may range from 10 km 
(6.2 mi) to perhaps 50 km (31 mi), and 

4. The large or synoptic scale which ranges from 50 to 
10,000 km (31 to 6214 mi). 

Winds on the storm or larger scales occur over such a long period, 
with respect to an aircraft's motion, that they are easily accommodated 
by the pilot and are of no concern other than how they might affect 
scheduled arrivals, fuel economy, etc. The turbulent scale accounts for 
bumpiness during aircraft flight, and is serious only when the "bumps" 
are intense, possibly resulting in structural damage to the aircraft or 
aircraft failure, excessive pilot work load, and passenger discomfort or 
injury. The scale distinction between turbulence and wind shear, 
however, is primarily a matter of interpretation. 

There is, however, substantial evidence that wind shears occurring 
on a scale of 1 to 10 km can create serious difficulties to aircraft, 
particularly in the landing or takeoff mode. The wind shear models 
discussed in the previous section incorporate scales of this magnitude. 
However, the high-frequency or high-wave number disturbances (i.e., 
turbulence) have been filtered out by the measurement technique or the 
extrapolation method inherent in the model. The wind speed components 
from the tower data are averaged over a 10-sec period. Thus, frequencies 
higher than 0.1 Hz are not contained within the measurement. For radar 
data, the wind speeds measured are the average value for a volume 
element typically 150 m long and of variable radial dimension. Finally, 
it is not clear what time or length scales are associated with wind 
speeds reconstructed from accident investigations; however, extrapo- 
lation of the flight data recorder values to a two-dimensional grid 
smooths the data immeasurably. 

To include the high-frequency components of thunderstorm wind 
variation (i.e., turbulence) current wind shear models generally super- 
impose simulated turbulence on the quasi -steady winds. The purpose of 
the turbulence is to insure a realistic pilot work load. The turbulence 
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simulation reported in Foy [19] is based on a Dryden spectrum with the 
intensity and length scale of each of the three velocity components 
programmed as a function of altitude. Unfortunately, the values of 
intensity and length scale utilized as input to these simulation models 
are not fully known for the thunderstorm or other severe wind shear 
environments. 

Measurements of the power spectral density function for turbulence 
in thunderstorms is reported as early as 1962 by Steiner and Rhyne [33]. 
Their data were measured over an approximate range of reduced spatial 
frequency of 0.004 to 0.4 rad/m. The theoretical von Karman spectrum 
follows the data in this frequency range very well as demonstrated by 
Houbolt et al . [34], see Figure 2.8a taken from this reference. The 
Dryden spectrum, on the other hand, does not compare as well with the 
data. All reported data were measured in the altitude range of 12 to 8 
km and thus are probably not representative of the low-level approach 
and takeoff environment. 

Houbolt et al . [34] gives a comparison of the power spectral 
density function of severe storms with that of cumulus clouds and clear 
air turbulence as shown in Figure 2.8b. One can see from Figure 2.8 
that the turbulence spectra for severe storms behave very similar to 
that of cumulus clouds and clear air turbulence with the only major 
difference being higher amplitudes of the power spectrum, which indi- 
cates higher turbulence intensity. 

Houbolt et al. [34] recommends for evaluation of the Dryden spec- 
trum a value of L = 1036 m (3400 ft) and values of a of 10.2 to 4.75 m/s 

(33.5 to 15.6 ft/s) for the vertical fluctuation and of 9.82 to 5.63 m/s 

(32.2 to 18.5 ft/s) for the lateral fluctuations. Again, these data are 

measured at very high altitudes and probably do not include effects due 
to the presence of the ground. The ground is expected, however, to have 
a strong effect on the turbulence length scale and intensities. No 
actual data for L and a below 500 m (4500 ft) nor how they vary with 
height in thunderstorm conditions is presently available. 

It should be noted that no turbulence information relative to 
the distribution of gust across an airfoil is available. These gusts 
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Reduced Frequency, rad/m 


(a) Measured and fitted^spectra 
for thunderstorm. L is 
integral length scale. 



(b) Typical power spectra of vertical component 
of turbulence measured in clear air, 
cumulus cloud, and thunderstorm. 


Figure 2.8 Turbulence power spectra for thunderstorm conditions [34]. 



can significantly influence the rolling motion of the aircraft and cause 
control upset to the point where corrective action may cause structural 
damage to the aircraft. A detailed discussion of turbulence modeling 
for thunderstorm wind shear is provided in Frost et al . [14]. 

2.4 Conclusions Relative to Wind Shear Models 

It is apparent that further data on wind shear is needed before 
standard models can be adopted for system qualification. Experiments 
are, however, in progress on gathering three-dimensional data for all 
three wind velocity components (Fujita et al . [32]) and for measuring 
the length scale and intensity of turbulence as well as distribution of 
turbulence across the airfoils during approach and takeoff in and near 
thunderstorms {Camp et al . [35]). Until these data are available and 
processed into the appropriate format for use in computer simulation and 
manned-flight simulator programming, the existing models are all that 
are available for use in qualitative analyses of aircraft performance in 
severe wind shear conditions. Their limitations, however, must be borne 
in mind. The following section describes the influence of wind shear on 
aircraft performance based on our existing knowledge of prevailing wind 
fields as described in the foregoing as well as on hypothetical models 
developed to isolate certain physical characteristics of wind shear. 
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3.0 AIRCRAFT PERFORMANCE IN WIND SHEAR 


3.1 Basic Considerations 

To investigate the effects of wind shear on aircraft performance, 
we must first examine the equations of motion of an aircraft in a 
variable wind field. The general form of the dynamic equations are, 
therefore, summarized in Appendix A. In the derivation of these 
equations, the earth is assumed to be a stationary plane in inertial 
space. This assumption is well justified for takeoff and landing 
problems, which are the main considerations in this study. Also, the 
airplane is treated as a rigid body having a plane of symmetry. This 
assumption implies that the motions of the atmosphere are of suffi- 
ciently large scale that they act uniformly over the airplane at any 
given moment. As noted earlier, little work has been done on the 
effects of wind shear distributed over the aircraft. This topic is 
addressed by Houbolt [36]. The spatial distribution of wind over 
aircraft has, however, been diagnosed as a significant factor in several 
recent aircraft accidents. For example, the National Transportation 
Safety Board (NTSB) reports that on February 24, 1980, a Beechcraft 
Bonanza BE-35 aircraft crashed near Valdosta, Georgia, during an encoun- 
ter with severe thunderstorms [37]. All the occupants were killed when 
their aircraft experienced an in-flight breakup. On August 26, 1978, 
two persons were killed when a Piper PA-28 aircraft experienced an 
in-flight breakup during an encounter with severe thunderstorms near 
Boulton, North Carolina. 

Also, large aircraft have experienced breakup during flight near or 
through thunderstorms [38]. A C-141 jet military cargo transport lost 
a wing when it broke up in flight over England. Although it is unlikely 
that the wing failure was caused by weather-induced stresses alone, 
reconstruction of the thunderstorm suggests very sharp downdraft gra- 
dients were encountered (wind speeds in excess of 51 m/s (100 kts) 
were suspected). 19 



Thus, it is apparent that strong wind gradients can pose a hazard 
to the structural integrity of an aircraft, as well as to its flying 
qualities. However, knowledge of the magnitudes of these wind shears in 
three dimensions is not yet available, and no analysis has been carried 
out to date on how such severe gradients would influence the lateral, 
roll, and yaw motions of the aircraft on approach and takeoff. To date, 
computer analyses and flight simulator studies are based on two- 
dimensional models of wind shear. Therefore in the following discussion, 
the major emphasis is placed on three-degrees-of-freedom analyses. As 
mentioned, studies are under way to measure the three-dimensionality of 
wind shear [32] or gust gradients across airfoils [39], and the results 
from these studies will be used to develop meaningful three-dimensional 
models of severe wind shear phenomena. 

As previously discussed, however, wind shear is defined in terms of 
relatively long-scale motions in the atmosphere. The higher frequency 
wind fluctuations (i.e., turbulence) may increase the pilot's work load 
but not necessarily affect the general flight path of the aircraft. 
Therefore, analyses based on the assumption that the wind fluctuation 
acts over the entire aircraft are expected to lead to meaningful quali- 
tative conclusions. 


3.2 Equations of Motion with Three-Degrees -of-Freedom 


To reduce the governing equations to three degrees of freedom and 
thus simplify the discussion as to how wind shear terms enter, consider 
Equation A. 2 in Appendix A. This equation reduces to the form: 

m(V + lI5 +qW )=T -D-mgsin0,, 

^ xw ^w ZVi' xw ^ w 


m[Q 


zw 


- q (V 


= ^zw 


- L + mg cos 0 


w 


(3.1) 


The nomenclature is given in Appendix B. From Equation A. 4a the wind 
components relative to the aircraft are related to those measured rela- 
tive to the earth by: 
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(3.2) 


“xw = “xE - “zE 


“zw “xE si" 8« * “zE "°s e^ 


Hence the time derivative becomes: 


'^xw = “xE ®w - '^zE si" s„ - (WxE si" ®w ^ “zE "°s 6„)e„ 
^zw = *^xE si" e„ + ^zE ®w * ^“xE "°s \ sin e„)e„ 


(3.3) 


where q = e . Thus, Equation 3.1 becomes 
w 


"'i' " i^xw ■ ii ■ ^i" ®w ■ '"^'^xE “S 6^ - sin e^) 

'"''®w " ‘"^zw + L - "la cos 0« '"^'^xE ^i" ®w '^zE ®w^ 


(3.4) 


Equation 3.4 shows the direct influence of wind shear terms on the rate 
of change of airspeed, V, and pitch angle, 9^. Of course, wind varia- 
tions also influence the values of D and L as discussed in Appendix A. 
Wind shear terms, however, do not appear directly in the moment equation 
(Equation A. 5), but wind shear does enter through the aerodynamic moment 
coefficients . 

Equation 3.4 expresses the influence of wind shear on the rate of 
change of airspeed and pitch referenced to the wind axes. These are the 
changes which the pilot observes from his airspeed indicator and flight 
director. The force equations can also be written in terms of the 
inertial velocity as 


= 

-L sin 6 

- D 

cos 6 

- mg 

sin 

®E "^XE 

mv^; = 

L cos 6 

- D 

sin 6 

- mg 

cos 

®E * ^zE 


In this case, wind shear terms do not appear directly in the equations; 
however, changes in wind are reflected not only in L and D but also in 
6, which is given by the expression 
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(3.6) 


6 = sin"^[(W^^ sin cos e^)/V] 

Regardless of whether Equation 3.4 or Equation 3.5 is used to 
compute the flight trajectory, the angle of attack, a, is directly 
influenced by wind shear through the relationship 

“w " ^ ^"^zw - L + mg cos sin cos e^)]/mV (3.7) 

This equation is obtained by rearranging Equation A. 7a. Details of the 
derivation of Equation 3.7 are given in Reference 40. 

3.3 Effects of Wind Shear Terms 


The effects of wind shear can be estimated by comparing the magni- 
tude of the terms appearing in Equation 3.4. Consider as an example an 
airplane having the characteristics of a B-727 descending at a sink rate 
of z = 3.75 m/s (7.3 kts) and a ground speed of x = 75 m/s (145 kts) 
in still air, the thrust (considered acting only along the x axis), 
lift, and drag per unit mass are approximately 

T/m = 1.14 m/s^(2.2 kts/s); 


L/m = 9.81 m/s'^(19.1 kts/s); 

D/m = 1.21 m/s^(2.4 kts/s). 

The derivative of the x-component of wind velocity for a wind 
varying only spatially is given by 

= x(8W^/3x)^ + z(9W^/9z)^ 

E 

and for the z-component by 


IaI = x(9W^/9x)r + z(9W^/9z)c 

Zr Z t Z t 


(3.8) 


Equation 3.8 is derived from Equation 3.3 where the subscript E is now 
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dropped for convenience and the approximation cos 0^=1, sin = 0, 

and 0, , = 0 is made, 
w 

Now consider Equation 3.1 with the same approximations: 


I _ D 
m ~ m 
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(3.9) 


For the terms on the right to be comparable in magnitude to T/m and D/m, 
3W /9x ~ 0.02 s"^ (1 kt/100 ft) and 3W /3z - 0.32 s“^ (19.5 kts/100 ft). 

X A 

In turn, for the wind shear terms in Equation 3.9 to be comparable in 
magnitude to L/m, then 3W^/3x - 0.13 s~^ (7.7 kts/100 ft) and 3W^/3z 
2.61 s"^ (159.1 kts/100 ft). 

Consideration of this simple calculation reveals that relatively 
large vertical gradients in the horizontal and/or vertical wind velocity 
components can be tolerated because the sink rate, z, during most 
approaches and takeoffs is small. Of significantly more interest is the 
observation from the simple calculation that relatively small values of 
shear in the horizontal direction result in values of the wind shear 
terms having the same order of magnitude as the lift and drag terms. It 
is obviously not the magnitude of the shear alone but the product of the 
horizontal velocity and the shear as well as the value of the glide 
slope, 0^, which dictates the strength of the wind shear effects. 

In practically all literature related to wind shear prior to 1977, 
magnitudes of vertical wind shear are reported. Values on the order of 
0.13 to 0.16 s"^ (8 to 10 kts/100 ft) are considered to be severe. 

These values correspond relatively close to the 0.32 s (19.5 kts/100 
ft) predicted by the simple calculation. It appears, however, that 
considerably more attention should be given to horizontal wind shear. A 
value of 8 kts/100 ft is the value of wind shear to which automatic 
landing systems are presently certified in the United States and the 
United Kingdom [41,42]. No discussion is given as to whether this is a 
vertical or horizontal shear. The above results suggest that 
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certification advisories should specify the value of wind shear which 
applies along the flight path. 

In the atmospheric boundary layer under normal conditions the major 
wind shear is in the vertical direction. The CAeM Extraordinary Session, 
1974, as reported in Reference 41, confirmed that statistics show there 
is approximately a 100 percent probability that the value of 8 kts/lOO 
ft will be exceeded on at least one landing per lifetime of the average 

aircraft. Reference 41 reports the following frequency for vertical 

wind shears of the given intensity: 

0.05 s"^ (3 kts/lOO ft) on 50 percent of the occasions 

0.08 s"^ (5 kts/lOO ft) on 17 percent of the occasions 

0.13 s”^ (8 kts/lOO ft) on 2 percent of the occasions 

0.16 s ^ (10 kts/lOO ft) on 0.4 percent of the occasions. 

There is, however, very little if any information relative to the 
expected frequency or intensity of wind shears in the horizontal direc- 
tion. This is partly due to the difficulty associated with measuring 
horizontal shears. Also, wind shears in the horizontal direction will 
be strongly influenced by terrain features, discontinuities in surface 
texture, and other microscale features. Fichtl et al . [43] summarized a 
number of surface features which can influence wind fields around 
airports. Frost and Camp [23] have also surveyed various meteorological 
phenomena which can create strong wind shear. 

From results reported by Goff [17], the average value of horizontal 
shears measured over several levels of a 500-m (1500 ft) tower during 
the passage of approximately 20 thunderstorms is 3W^/8x = 0.09 s“^ (5.3 
kts/lOO ft) and 8W /9z = 0.04 s“^ (2.4 kts/lOO ft). These values were 

A 

computed with Taylor's hypothesis, which implies that the variations in 
wind are frozen in the flow field as the storm passes over the tower. 
Hence, the values are considerably smaller than instantaneous values. 

It is readily apparent, however, that horizontal wind shears in thunder- 
storm conditions very quickly exceed the magnitude of 0.02 s"^ (1 kt/100 
ft) estimated as significant from simple calculations. 


24 



Table 3.1 shows the lift and drag to mass ratios for three aircraft 
types. These values are for roughly steady flight at landing speeds in 
which case the lift per unit mass is very close to the value of the 
acceleration of gravity, g. Thrust or drag per unit mass depend upon 
the lift/drag ratio. If this ratio is high, then thrust is low, and the 
effect of the terms x8W /9x and z9W /9z — which act like horizontal 

A A 

forces — will be relatively large (see Equation 3.1), This is all the 
more true if the aircraft is travelling at high speed. The above argu- 
ments suggest commercial airliners which have lower values of D/m and 
higher values of are more susceptible to wind shear than smaller, 
lighter aircraft. 


3.4 Qualitative Analysis 

The influence of wind shear on aircraft motion is described quali- 
tatively in several recent articles [44-53]. Melvin [44] appears to 
have given one of the first descriptions of the effect of wind shear 
during approach; this is summarized as follows. 

When a wind shear is encountered during approach, the effects are 
twofold and opposite in direction. One effect is dependent upon the 
rate of the shear while the other is dependent only upon the magnitude 
of the shear. 

The first effect is associated with the attempt to maintain a 
prescribed airspeed. If an aircraft is on an approach at 62 m/s (120 
kts) IAS with a 10.3 m/s (20 kt) head wind, ground speed will be 
51 m/s (100 kts). If the head wind ceases, the aircraft will need to 


TABLE 3.1. Typical Values of V^, D/m, and L/m of Different Aircraft Types 


V^, m/s (kts) 
D/m, N/kg 
L/m, N/kg 


DHC-6 

46 (90) 

3.420 

9.807 


B-727 

72 (140) 
1 .008 
9,807 


Queen Air 

56 (110) 

1.230 

9.807 
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accelerate to a ground speed of 62 m/s (120 kts) to maintain its 
airspeed. This can be accomplished by pushing the nose over and accept- 
ing a loss of altitude or by prompt application of thrust to accelerate 
the aircraft at a rate equivalent to the rate of wind shear. 

The second effect is associated with the attempt to fly a pre- 
scribed glide slope. Consider an aircraft flying a 3° ILS on a stabi- 
lized approach. If the aircraft described above encounters instan- 
taneous wind shear from a 10.3 m/s (20 kt) head wind to no wind, the 
airspeed will drop from 62 m/s (120 kts) to 52 m/s (100 kts), the 
nose will pitch down, and the aircraft will drop below the glide slope. 
The loss in altitude will be directly proportional to the new wind 
condition, assuming thrust is maintained constant (the principle of 
exchange of potential energy for kinetic energy). Once the energy 
exchange is accomplished, the aircraft has more thrust than is required 
to fly the glide slope under the no-wind condition. Thus, it will 
gradually gain on the glide slope and overfly it. 

The apparent effect of a decreasing head wind is illustrated in 
Figure 3.1. The effect is different depending upon where the shear 
occurs relative to the ground, the rate of shear, and the magnitude of 
shear. If the wind shear occurs very close to the ground, the aircraft 
will hit short. On the other hand, if the shear occurs some distance 
from the ground, the aircraft will tend to overfly the touchdown zone. 
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Figure 3.1 Head wind shearing to tail wind or calm [48]. 
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As Melvin points out, however, no pilot should attempt to maintain 
glide slope with a constant thrust setting [47]. In the high head wind 
before encountering the shear, the pilot will be using a larger thrust 
setting than is required to fly the glide path in a no-wind condition. 
When the wind begins to decrease, the aircraft will tend to lose air- 
speed and fall below the glide slope. The pilot, in recognizing this, 
will add thrust to return to the glide path. (Theoretically, the amount 
of thrust required to equal that required to accelerate the aircraft 
mass at the same rate the wind is shearing.) Once the aircraft is back 
on the glide slope, the pilot will need to gradually reduce thrust to 
account for the lessening head wind. When the wind shear ceases, the 
aircraft no longer needs to accelerate and a thrust reduction should be 
applied to prevent overflying the glide slope. 

Figure 3.2 illustrates the condition of a tail wind which rapidly 
decreases to a calm or head wind. Initially, the IAS and pitch will 
increase and the aircraft will overfly the glide slope. To compensate 
for this, a thrust reduction is required initially to reduce the air- 
craft's high ground speed, followed by a gradual thrust increase. When 
the wind ceases altogether or changes to a head wind, a large thrust 
addition is required to restabilize power after the initial reduction 
and to prevent loss of ground speed. 

Once again, the effect of the tail wind shearing to a calm or head 
wind is dependent upon the altitude at which the shear occurs. If the 
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Figure 3.2 Tail wind shearing to head wind or calm [48]. 
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shear occurs close to the ground and the thrust is not reduced quickly 
enough, the approach will be high and fast with the danger of over- 
shooting. If, on the other hand, the shear occurs well above the 
ground, the aircraft will first rise above the glide slope and, if the 
thrust is held relatively constant, sink back below the glide slope, 
landing short. 

Higgins and Patterson [51] have also looked qualitatively at flying 
procedures in hazardous wind shear. They used static performance curves 
to provide pilots with some ideas relative to handling shears. They 
point out that if implemented, these ideas would aid in avoiding catas- 
trophe if the pilot's aircraft was inadvertently caught in a combination 
of severe downdraft and/or severe wind shear that resulted in high rates 
of descent and/or severe loss of airspeed, especially within approxi- 
mately 122 m (400 ft) of the ground. They also discuss the following 
points: 

t Basic performance conditions 

• Airplane energy management concepts 

f Maneuvering margins 

f Angle of attack consideration 

• Attitude considerations 

t Performance effects of acceleration along the flight path 

• Performance effects during flap retraction. 

The key points of some operational techniques they recommend relative to 
hazards of landing, approach, and takeoff in wind shear environments 
are: 

• When forced to fly at speeds near stick-shaker because of 
wind shears, good climb performance and maneuver margins 
still exist. Rapidly accelerating the aircraft away from 
stick-shaker could result in a significant loss of altitude. 

• High attitudes are required at stick-shaker speeds and go- 
around thrust to attain the maximum climb capability of 
the aircraft. 
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t Rapidly accelerating to maintain Vref V 2 airspeeds 
during a wind shear will severely reduce climb capabili- 
ties. Conversely, decelerating to stick-shaker speeds can 
provide added climb capability to compensate for large 
downdrafts. 

The recommendations are based on performance analyses from charts 
which are valid for stabilized 1-g flight conditions at constant indi- 
cated airspeeds for the airplane weight given in the report. The 
authors point out that if pilots make use of any of the specific atti- 
tudes from these charts as a guide for operation of a B-727, the 
attitudes should be treated only as initial targets. Flight in severe 
wind shear is a dynamic, constantly changing situations and confirma- 
tion that any given attitude is adequate for any given situation comes 
from instrument readings which show that the aircraft is responding in a 
satisfactory and desirable manner. 

ALPA's Airworthiness and Performance Committee (see Steenblik 
[52]), on the other hand, is concerned that many airline flight training 
departments continue to train pilots to promptly trade airspeed for 
altitude by pitching up until the airspeed decays enough to activate the 
stick-shaker (last recommendaton by Boeing article). The committee 
argues that pilots should attempt to achieve minimum drag speed (best 
angle of climb speed) during wind shear encounters. 

When performance is critically limited by wind shear effects, the 
ALPA committee recommends that pilots fly near the minimum drag point 
for best climb angle performance until ground impact is eminent, then 
exchange all available energy to flair the aircraft and soften the 
impact or to sustain flight in ground effects until clear of the shear. 

A distinction must be made between excess thrust over drag capa- 
bility, which contributes to long-term flight path performance and 
energy trades (kinetic for potential and vice versa). A turbojet 
aircraft attains its maximum climb angle performance at approximately 
its minimum drag speed. There is a small range below the exact minimum 
drag speed for which drag does not increase significantly, but drag does 
increase rapidly as speed is lowered, rapidly reducing climb angle 
performance. 
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It is unreasonable to think that any pilot would deliberately fly 
up the back side of a drag curve (see Figure 3.3) when performance is 
limited by wind shear. However, this does not mean that a pilot, upon 
realizing that impact with the ground or an obstacle is eminent, would 
not pull the aircraft up and sacrafice airspeed to avoid or reduce 
impact. The best climb performance occurs when the aircraft is most 
energy efficient and that is at the minimum drag point. Aircraft per- 
formance limited by wind shear cannot be increased by flying up the back 
side of the drag curve. 

Steenblik [52] concludes that in the 1975 Continental Airlines' 
takeoff accident in Denver, there was no way the aircraft could acceler- 
ate inertial ly fast enough to overcome the effects of the shear and 
avoid ground contact. Steenblik [52] believes that it would have been 
impossible for the pilot to have flown out of the shear conditions by 
just increasing the pitch. The ALPA committee goes on to point out that 



Airspeed (kts) 

Figure 3.3 Flaps 30 constant speed climb capability [51]. 


30 



there appears to be some confusion relative to the all -engine climb 
capability of an aircraft at the stick-shaker speed. Considering that 
the climb capability might be approximately 1200 ft/min for a represen- 
tative aircraft, many fail to recognize that the all-engine climb 
capability is probably double that amount at the minimum drag point. 
Deliberately trading all available energy down to the stick-shaker speed 
while increasing drag to the point of drastically reducing climb capa- 
bility is, in the committee's belief, an unsafe practice. 

If an aircraft is operated at the minimum drag point (or in the 
fairly flat portion of the drag curve), it will achieve its best perfor- 
mance; then if energy is traded as the aircraft enters ground effects, 
the trade could result in a successful go-around, but at least would 
result in minimizing ground impact. If the initial climb moves the 
aircraft well up the back side of the drag curve, however, there will be 
reduced capability for a sustained climb and energy available as a last 
resort will be very limited. 

The importance of carrying extra speed for landing cannot be 
overemphasized; such an energy trade can give much faster results than 
improved performance from an increase in thrust. An energy trade, 
however, is a one-shot affair; only a thrust increase can make a long- 
term contribution to a new flight path. Pilots should not be reluctant 
to trade energy down to the minimum drag point — after that point, energy 
to be traded should be reserved to use when ground impact is eminent. 

Excess energy to trade down to the minimum drag point is important, 
as is the early recognition of the effects of wind shear and the rapid 
application of maximum thrust. Bliss [53], however, concludes that the 
airspeed/ground speed concept is essential for wind shear protection. 

He believes that without the correct ground speed value, a wind shear 
warning system cannot solve the problem adequately, and no amount of 
training can be of any use for a severe wind shear situation. This, he 
claims, is true for either a manual or coupled approach. The accelera- 
tion margin device, in Bliss's opinion, also becomes useless on a non- 
precision approach. 
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A pilot's conception of speed is traditionally oriented solely with 
reference to indicated airspeed (IAS). The conventional performance 
charts (those used by Higgins and Patterson [51]) are therefore refer- 
enced to indicate airspeed values which Bliss believes are valid only in 
static air, but worthless in a wind shear. In every case where the 
control of an airplane is placed in a hazardous condition due to wind 
shear, it is specifically the result of either an excessively high or 
excessively low ground speed value at an attitude low enough to compro- 
mise or preclude recovery. 

To establish the magnitude of the effects described above, analyt- 
ical models with varying degrees of sophistication have been developed. 

The author is unaware of any models, however, which include all six 
degrees of motion with any realistic three-dimensional model of severe 
wind shear. Thus, only results for three-degrees-of-freedom and two- 
dimensional wind fields are discussed in the following section. 

3.5 Mathematical Analysis 

Etkin [54], in 1946, appears to be one of the first to analyze 

flight in wind shear. Using a system of linearized equations, he 

investigated the performance of a light airplane gliding at 27 m/s (52 

kts) through wind shears of 3W /3z = 0,04, 0.03, and 0.002 s”^ (2.6, 

/\ 

1.7, and 0.9 kts/100 ft). He predicted that the aircraft would over- 
shoot the desired touchdown point by roughly 792, 549, and 274 m (2600; 
1800, and 900 ft), respectively. These results are based on the assump- 
tion of an approach at constant relative velocity, V, and constant pitch 
angle, y. The results of Etkin 's study for the magnitude of overshoot. 

Ax, can be expressed as 

AX = rz?/2Vg sin y^ (3.10) 

where r = 3W^/3z and Ax = is the length of overshoot. 

Etkin concludes that the distortion of the flight path during both 
glide and climb is greatest when the rate of descent is small and when 
the wind velocity is large relative to airplane speed. Thus, the 
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aircraft: most affected would be machines with low wing loading and flat 
glides. Fast machines with steep glides would be less disturbed. This 
result appears to be in variance with the conclusions reached above. It 
is shown later that flight path calculations using a linear wind shear 
in the vertical direction give seemingly contradictory results to the 
same computation using a logarithmic wind speed profile, which is more 
realistic of the atmospheric boundary layers. 

Gera [55], using a similar system of equations as Etkin, assessed 
the influence of wind shear on the longitudinal motion of airplanes. He 
also assumes completely horizontal wind with linear variation in speed 
with altitude and arrives at almost the identical conclusions as Etkin. 

Gera shows that the deviation from the touchdown point in wind 
shear as contrasted with steady winds is expressed by the relationship 


Ax 


g 

Sin Ay^ 



sin Ay^ 
sin Tg 


(3.11) 


where a is a nondimensional wind shear parameter defined by 
O = Vg(3W^/3z)/g (3.12) 


This expression gives the amount of overshoot and undershoot at ground 
level relative to the flight path and steady wind. Gera concludes that 
in a head wind decreasing with altitude there is an undershoot (ax < 0) 
as long as the inequality 




2V? 

e 

sin Ay 

^ 

sin Yg 


(3.13) 


is true. If the height lost during the descent through the shear layer 
exceeds the right-hand side of the inequality, an overshoot will occur. 
It is important to note, however, that for moderate values of wind shear 
there is always an initial undershoot regardless of the thickness of the 
shear layer. This result is in agreement with the qualitative discus- 
sion given earlier. 
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In Gera's analysis for the above result, the undershoot was calcu- 
lated for the case in which the airspeed, angle of attack, and throttle 
setting were the same in both steady wind and wind shear. It was found 
that these conditions were possible if the airplane in wind shear 
assumed a pitch attitude different from the steady wind value. For the 
case where the airspeed, angle of attack, and pitch angle are the same 
for both the conditions of steady wind and wind shear (possible, at 
least in theory, by having different control deflections and throttle 
settings under the two conditions). Equation 3.11 reduces to Equation 
3.10. The value of Ax in this situation for normal glide is always 
positive, and hence an overshoot always occurs. The reference condi- 
tions for Gera's analysis are zero external moment and constant air- 
speed, angle of attack, and pitch angle. Etkin began his analysis with 
an arbitrary initial value of angle of attack. This value became 
constant, however, in a very short period of time; so essentially, the 
analysis was for constant angle of attack (which in this case is also 
constant pitch angle) and constant airspeed. 

Gera [55] also investigates the effect of linear wind shear on 
longitudinal stability. For the assumption of constant airspeed, 
constant angle of attack, and pitch angle, the effect of wind shear on 
the short-period motion and the phugoid damping was negligible, but the 
phugoid frequency and damping ratio were found to vary considerably with 
wind shear. The time for the phugoid to damp to half amplitude increases 
in a climb and decreases in a dive, as expected. Positive shear (head 
wind changing to a tail wind, which is typically the situation associated 
with flying through a downburst) is shown to amplify these effects. He 
also notes that wind shear affects the phugoid mode even in level 
flight. 

Etkin [4] analyzes the longitudinal stability of a typical STOL 

“5 ”1 

airplane for linear vertical wind shears from -3.4 x 10 s (-0.002 
kts/100 ft) (the head wind case) to 3.4 x 10 ^ s ^ (0.002 kts/100 ft) 

(the tail wind case). These are very low wind shear values. He found, 
however, the effects on both the phugoid and pitching mode to be very 
large. A strong head wind decreases both the frequency and damping of 
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the phugoid, and a strong tail wind changes the real pair of matching 
roots into a complex pair representing a pitching oscillation of long- 
period and heavy damping. 

Sherman [56] also used a linearized system of equations to carry 
out a stability analysis of wind shear effects on a large jet transport. 
In the case of the phugoid mode, positive wind shear (a shear that 
changes head wind to tail wind) caused the phugoid to remain periodic 
and stable, although the time to damp to half amplitude decreased as the 
shear gradient increased. A negative shear (a shear that changes tail 
wind to head wind) caused the phugoid to become unstable for values of 
the shear parameter, a, greater than unity (i.e., a > 1). 

The general conclusions from the above are that wind shear has a 
pronounced effect on the phugoid modes of aircraft stability but little 
or no effect on the short-period modes. The paper by Moorhouse [57] 
lends further support for this argument. 

The assumption of linear wind shear discussed in the preceding 
results does not present a realistic simulation of the atmospheric 
boundary layer wind profiles. Typical wind speed profiles found in the 
atmospheric boundary layer under moderate climatological conditions are 
best represented by a logarithmic profile [58,59] 

Wx = ^ [nn(z^/Zg) + .(.(z^/L)] (3.14) 

The function i|j(z^/L) (note L is the Monin-Obukhov length scale) takes 
different forms depending on the characteristics of the atmospheric 
boundary layer, i.e.. 

Neutral Boundary Layer : 

^{z^/D = 0 (3.15) 


Unstable Boundary Layer : 


z^/L 


-1/4- 


ip(z^/L) = z^/L[l - (1 - 18z^/L)-''^]d(z^/L) 




( 3 . 16 ) 
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(3.17) 


Stable Boundary Layer : 
i|i(z^/L) = 5.2z^/L 

These forms of wind profiles do not lend themselves to linearized 
models. 

Luers and Reeves [60] attack this problem by utilizing a nonlinear 
system of equations similar to Equation 3.5. The landing of seven com- 
mercial/mil itary-type aircraft was computer simulated starting from an 
initial altitude of 90 m (300 ft). Deviations in touchdown point in 
excess of 910 m (3000 ft) resulting from variation of the horizontal 
wind during the final 90 m (300 ft) of descent with no pilot or auto- 
pilot feedback were computed. Vertical wind shear associated with a 
range of values of surface roughness, z^, and stability, L, were 
considered. 

Their analysis, however, neglected the wind shear terms in the 
rate-of-change-of-angle-of-attack relationship given by Equation 3.7. 
With these terms included, it is found [61] that the results of Luers 
and Reeves overpredict the magnitude of the touchdown deviations by 
roughly 60 percent. The trends of their results are correct, however, 
even though the magnitude of deviation from the desired touchdown point 
is too high. 

Interestingly, the results of the flight path analysis utilizing 
the nonlinear system of equations and the logarithmic wind shear pro- 
files resulted in the aircraft undershooting the runway when landing in 
a head wind and overshooting the runway when landing in a tail wind. 

This result is in direct contrast to the results reported by Etkin [54] 
and Gera [55] for linear wind shear profiles. Undershooting the touch- 
down point was also reported by Frost and Reddy [62] and by Denaro [63]. 

Denaro [63], analyzing aircraft flare, explains this apparent 
contradiction in the effect of wind shear by a combination of two 
factors. First, a logarithmic wind shear has a higher rate of change of 
wind magnitude at the lower altitudes than does a comparable linear 
shear. (Note the wind shear for a logarithmic profile goes as 1/z^ and 
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becomes very large near the ground.) An aircraft at flare altitude is 
likely to experience large wind gradients and therefore large airspeed 
changes. In a logarithmic head wind shear, the aircraft is signifi- 
cantly below the normal speed for which the flare control is designed. 
Second, because the onset of shear is gradual, the throttles do not 
respond as much initially as they do in the linear shear case. Conse- 
quently, when the aircraft reaches flare altitude in a logarithmic head 
wind shear, it does not have greatly advanced throttles as in the linear 
head wind shear case. Therefore, throttle freeze and retard at this 
point have a significant effect on reducing airspeed even further. The 
aircraft starts flare at a rather low sink rate, but response is poor 
and the sink rate is not well arrested. With the higher sink rate in 
the latter stages of flare, the aircraft lands both hard and short. 

Frost and Reddy [62] and Luers and Reeves [60], however, found 
short landings even without a control system being involved. In this 
case the difference between the results for a linear profile and a 
logarithmic velocity profile are a consequence of the initial trim 
conditions used to start the computation. For a fixed control system, 
the aircraft is trimmed at the value of wind shear existing at the 
initial altitude from which the calculation begins. With a logarithmic 
velocity profile at sufficiently high altitudes, the wind shear is very 
low. As the aircraft approaches the ground, the wind shear for a 
logarithmic velocity profile increases rapidly whereas the linear wind 
shear remains constant. However, the thrust and elevator setting have 
been set for the lower magnitude wind shear. Inspection of Equation 3.4 
illustrates that for fixed thrust and possible decreasing values of drag 
and lift due to reduced wind speed (note that L and D are functions of 
angle of attack and the rate of angle of attack as well as other tran- 
sient variables [40]), the increasing wind shear term strongly influences 
the sink rate of the aircraft. Frost and Reddy [62] had no difficulty 
removing this fast sink rate when an automatic control system was 
incorporated into the computer analysis. 

The preceding analyses have investigated only shear of the hori- 
zontal wind. Under more severe wind shear conditions, particularly 
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thunderstorms, major fronts, and flow fields around buildings or other 
surface terrain features, the vertical wind component can be extreme. 
References 40, 62, 64, 65, 66, and 67 report investigations with the 
vertical wind speed component included in the equations of motion. The 
impetus to investigate flight with severe variations in both vertical 
and horizontal wind speeds was generated by the Eastern 66 accident in a 
severe thunderstorm at JFK International Airport on June 24, 1975. This 
accident created immense concern relative to flight through thunderstorms. 
Frost et al. [16] and Foy [19] developed mathematical models of wind 
fields associated with strong environmental shears. Frost and Crosby 
[40] and TurkeT and Frost [64], utilizing these models in the form of 
computer table lookup routines, investigated the flight of various types 
of aircraft through thunderstorms and other strong wind shear condi- 
tions. The following sections describe these results along with results 
from other studies reported in the open literature. 
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4.0 FLIGHT IN STRONG WIND SHEAR ENVIRONMENTS 


4.1 Fixed Control Models 

Initial numerical studies of flight in thunderstorm-type wind shear 
were carried out under the assumption of fixed controls. Figure 4.1 
shows the computed descent of a DC-8-type aircraft through 11 different 
thunderstorms. In nearly all of these approaches an oscillation near 
the phugoid frequency of the aircraft is strongly amplified. This 
directly supports the conclusions of the stability analyses relative to 
the phugoid mode described earlier. 

McCarthy and Blick [66] independently analyzed flight in thunder- 
storms. Using a linearized model and a superposition technique, they 
investigated the flight behavior of a B-727-type aircraft in a thunder- 
storm, using wind data that had been obtained from in-flight measure- 
ments near thunderstorms. They also found amplified flight path 
oscillations at frequencies near the phugoid frequency of the aircraft. 

Frost and Crosby [40] applied their models to a number of aircraft 
types. The flight paths computed for two of the more severe thunder- 
storms are shown in Figure 4.2. In nearly all cases, the aircraft 
demonstrated high-amplitude oscillations at frequencies near the phugoid 
frequency. Table 4.1 shows the computed phugoid period and horizontal 
wavelength versus those predicted by simple theory. For the commercial - 
type aircraft, the frequency of the oscillations observed in the thunder- 
storm flight paths are very close to those predicted for the phugoid 
oscillations from simple theory. For the smaller DHC -6-type aircraft, 
the oscillations occurred at a somewhat higher frequency than the 
classical phugoid frequency. Correspondingly, the smaller aircraft 
showed less sensitivity to the thunderstorm wind fields. These results 
suggest that thunderstorm wind fields have characteristic scales of wind 
shear which can create hazardous oscillations in the flight paths of 
commercial -type aircraft. Severe oscillations in airspeed were also 
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Figure 4.1 Flight paths of DC-8-type aircraft landing with fixed 
controls at a -2.7° glide slope (numbers on curves 
designate different thunderstorm cases). 



z/h 

Figure 4.2 Comparison of different types of aircraft landing with 
fixed controls in thunderstorm cases 9 and 11 at a 
-2.7° glide slope. 
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TABLE 4.1. Phugold Period and Horizontal Wavelength. 



V 

T (sec) 

X 

(m) 

X = 

X/h 

Aircraft 

Type 

m s"^ 

Com- 

puted 

Pre- 

dicted 

Com- 

puted 

Pre- 

dicted 

Com- 

puted 

Pre- 

dicted 

DC -8 

70 

29.9 

31.7 

2,180 

2,203 

23.84 

24.09 

B-747 

66 

28.8 

30.0 

2,067 

2.085 

22.60 

22.80 

DHC-6 

46 

27.1 

20.7 

2,405 

1 .016 

26.3 

11.11 


T = /^V/g; X = VT 


computed. Figure 4.3. McCarthy et al . [67,68] computed almost identical 
results with their model. 

Augmentation of the phugoid mode during flight through severe wind 
shear suggests an accident-causing factor. McCarthy et al . [68] states 
that longitudinal wind gusts providing energy at the phugoid frequency 
may result in airspeed oscillation of a nature that would be difficult 
to control and, in fact, may lead to stalls and otherwise disastrous 
results. 

Most pilots are adamant that because of the low frequency of the 
phugoid oscillations, they can be controlled without difficulty. 



Figure 4.3 Comparison of indicated airspeed of DC-8-type and B-747- 
type aircraft landing with fixed controls in thunderstorm 
case 9. 
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However, during approach through a thunderstorm with other distractions 
such as poor visibility, runway slipperiness, etc., the effects of the 
phugoid oscillations can be insidious, and before the pilot realizes the 
presence of the oscillations, he may have reached a situation that is 
uncontrollable. In actual fact, this effect of first rising and then 
falling below the glide slope is exactly that described qualitatively by 
Melvin [44,45]. 

Additionally, nearly all conclusions relative to the phugoid oscil- 
lations by pilots and aerodynami cists alike are not based on the concept 
or experience of a forcing function (i.e., variable wind speed) driving 
the system at its critical frequency. In fact, most training is based 
on steady winds. In turn, there is some question as to how well phugoid 
oscillations are reproduced in a manned flight simulator. Thus, it is 
believed that the effect of forcing the aircraft at its phugoid fre- 
quency can be hazardous and should not be taken glibly. 

In many of the thunderstorm analyses carried out [26,69], the 
concept of a downburst, or extreme downdraft in the heart of the thun- 
derstorm cell, is suggested as the significant wind component contrib- 
uting to loss of flight control. Thus, the vertical wind is considered 
the prime factor creating hazardous conditions. Nearly all of the 
previous arguments, however, suggest that the horizontal wind component 
is equally important in creating flight hazards. To test the individual 
effects of the wind speed components, the computer program was run first 
with only the longitudinal wind component and then second with only the 
vertical wind component. Figure 4.4 shows the separate effects of the 
two wind components for a DC-8-type aircraft landing with fixed controls 
in a typical thunderstorm outflow. It is apparent that in the absence 
of the longitudinal wind component, the influence of the wind on the 
aircraft flight path, is considerably reduced. McCarthy et al. [68] 
arrived at identical results. 

4.2 Automatic Control Systems 

The preceding results show that serious departures from the glide 
slope occur during simulated landing of aircraft with fixed controls in 
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Figure 4.4 Comparison of DC-8-type aircraft landing with fixed 

controls in thunderstorm case 9, considering individual 
wind components separately and combined. 


thunderstorm gust fronts. Since the assumption of fixed controls is not 
realistic. Frost and Crosby [40] investigated automatically controlled 
flight. The automatic control systems using variable gains almost 
completely eliminate the severe perturbations from the flight path for 
the thunderstorm models considered in the study, Figure 4.5. However, 
the large control inputs and small response times required for the auto- 
matic control system to track the glide path in the thunderstorm cases 
may be difficult to achieve in operational hardware. 

Figure 4.6 shows the thrust control necessary to maintain the glide 
path during approach through a thunderstorm. The thunderstorm studied 
resulted in a tail wind shearing to a head wind. The insert in Figure 
4.6 shows the correspondence between Melvin's [44] qualitative descrip- 
tion of thrust requirement and that predicted by the computer simulation. 

4.3 Pilot Model s 

The automatic control computer model [40,62] was then expanded to 
incorporate a simulation of a human pilot into the computed response of 
the aircraft in wind shear [64]. 
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Figure 4.5 Flight path comparison of DC -8-type aircraft landing with 

(1) fixed controls, (2) automatic controls, and (3) automatic 
controls with turbulence included, in several different 
thunderstorm cases. 


Human pilot transfer function data were taken during compensatory 
tracking simulator experiments by Adams and Bergeron [70]. They tested 
six pilots (ages 30 through 47) and two test engineers in a flight 
simulator equipped with an oscilloscope and control stick. The sub- 
jects' static gains, lead and lag time constants throughout the runs 
were measured and variations between the subjects for given controlled 
dynamics (degree of vehicle controllabil ity) were examined. 

For the eight pilots the values of the transfer functions using a 
response time step of 0.01 second ranged between 0.254 and 0.905. For the 
reported study of Turkel and Frost [64] , a parametric study of pilot 
performance ratings between zero and one was considered more useful than 
using any specific pilot rating from the data of Adams and Bergeron 
[70]. The pilot model was incorporated into the control loop of Frost 
and Reddy [62]. Since the pilot does not move the servos to correct the 
plane's deviations as efficiently as the autopilot and, in effect, 
always lags the autopilot, the pilot's control signal inputs were 
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Figure 4.6 Rate of change of thrust required of DC-8-type aircraft 
landing with an automatic control system in thunderstorm 
cases 9 and 11 . 

reduced by a "perfection percentage," where 0 percent corresponds to 
zero inputs (fixed stick), 100 percent corresponds to "perfect" auto- 
pilot control, and a rating of 50 percent can be said to be average. 

This was accomplished in the program by multiplying the control signals 
by the "perfection percentage." The pilot constantly attempts to return 
the aircraft to the desired state but this occurs at a slower response 
rate than the "near-perfect" automatic control system. 

Fixed-stick, autopilot, and manned performance were compared for a 
B-727-type medium-sized commercial transport and for a Queen Air small 
commuter-type aircraft flown through a glide slope longitudinal wind 
profile detected by Doppler radar [71]. The wave-form wind disturbance 
was shown to excite the phugoid oscillations of both aircraft when they 
were flown in the fixed-stick mode, but presented no control problems 
for manned aircraft. 
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To investigate the significance of wind shear with a frequency 
equal to the phugoid, a fictitious quarter-sinusoidal wind field was 
modeled. A simulation was made for fixed-stick, autopilot-controlled, 
and manned aircraft with characteristics of a B-727 through this profile 
for a 6 m/s (12 kt) amplitude head-wind-to-tail -wind phugoid-frequency 
shear wave. This case revealed phugoid oscillations but clearly showed 
that this shear wave was not a serious problem for a manned vehicle. 

However, in a stronger disturbance — 10 m/s (19 kt) head-wind-to- 
tail-wind phugoid-frequency shear wave — significant deviation from the 
glide slope was noted for the autopilot, the 50-percent-rated pilot, and 
the 25-percent-rated pilot flight simulations, although no hazardous 
situations occurred. The low performance 5-percent-rated pilot initially 
lost control of the aircraft and dropped farthest below the glide slope. 
However, thrust was eventually increased to bring the aircraft back to 
the gl ide slope. 

In flight simulations through a full 14 m/s (27 kt) phugoidal- 
frequency sine wave, comparisons were made between autopilot control and 
control by pilots of varying skill. The autopiloted aircraft executed 
the best approach, while the high-skilled pilot descended below the 
glide slope but was eventually able to bring the aircraft back onto the 
glide path. However, the low-skilled pilot could not maintain adequate 
control and landed short. 

Classifying pilot response by means of a performance rating encom- 
passes the many intangibles encountered in pilot modeling which are too 
complex to simulate. These intangibles include pilot personality, 
training, knowledge, and warning of the encountered wind shear, as well 
as the element of surprise. Hence, a pilot with a low performance 
rating (for example, 0.03, which corresponds to a minimal control input) 
may be classified as poorly trained, slow-to-react , unknowledgable, or 
uninformed of the eminent wind shear. The report concludes that more 
work is clearly needed on pilot modeling, specifically to determine 
pilot response to the wind shear environment. However, for purposes of 
the study, the proposed pilot's "perfection percentage" gave useful 
results. 
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4.4 Comparison of Computer Simulation with Manned 
Flight Simulator Studies 


4.4.1 Description of Stud y 

Frost et al . [72] and McCarthy and Norviel [73] compared computer 
simulation with manned-flight simulator studies. The aim of this work 
was: 

1. To utilize the three-degrees-of-freedom aircraft trajec- 
tory computer program to examine aircraft/pilot response 
through wind shears including longitudinal sine waves, 

S-shaped waves, 1 - cosine vertical winds, and combina- 
tions at various frequencies and amplitudes as approxi- 
mations to the winds encountered in a thunderstorm 
downburst cell. 

2. To determine if the control system algorithm and aircraft 
trajectory program combination gives an accurate repre- 
sentation of the behavior of the real pilot by comparing 
the computed results with those measured in a manned 
flight simulator when subjected to the same input wind 
field models. 

The aircraft computer program was a three-degrees-of-freedom 
(horizontal, vertical, and pitch) program. This program and the pilot/ 
control system models are described in detail in Turkel et al . [74]. 

For comparison purposes, quantitative flight path deterioration param- 
eters were defined. These parameters were investigated to determine the 
degree to which they serve as a measure of hazardous flight conditions 
existing on approach through sinusoidally varying winds. The sinusoidal 
winds are an idealization of winds associated with flight through 
thunderstorm cells. 

4.4.2 Idealized Wind Speed Profiles 

The wind speed profiles selected for study are shown in Figure 4.7. 
They are based on the observation that an aircraft flying through a 
downburst would first encounter an increasing head wind with the wind 
changing to zero and resulting in an increasing and then decreasing tail 
wind. Depending upon the wind storm, this may either have a full sine 
wave effect or a S-shaped or half-sine-wave effect. In turn, the air- 
craft would encounter an increasing downdraft reaching a maximum at the 
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used to simulate a thunderstorm downburst cell. 


center of the downburst and then decreasing again to zero at the far 
side of the outflow. Evidence of this type of wind field has been 
determined by a number of studies (see Section 2.2). Figure 4.8 shows 
the average wind speed for 20 thunderstorm cases along a 3° glide slope 
at three different elevations. Note that the horizontal wind shear 
clearly illustrates an S-shaped sinusoid while the vertical wind demon- 
strates a similar profile to the 1 - cosine shape shown in Figure 4.7. 

Goff [15] examined the periodic nature of thunderstorm data. He 
computed the wind shear energy spectrum for the longitudinal wind compo- 
nent and the vertical wind component. Figure 4.9 shows that the energy 
for thunderstorm wind shears is contained in a frequency range that 
encompasses the typical phugoid frequency of most aircraft. The scale 
across the top of the figure indicates that the peak in the energy 
spectrum occurs somewhere near 100 to 50 seconds with respect to the 
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a) spectrum (longitudinal component) b) spectrum (vertical component) 

Figure 4.9 Wind spectra indicating frequencies associated with 

thunderstorm wind shear (n - number of observations in 
the time series and c = the mean tower wind speed) [15], 





aircraft's motion. This is equally true of the vertical spectrum, 
although note that the vertical spectrum contains very little energy 
when compared to the longitudinal component. Finally, qualitative 
inspection of Figure 2.6, page 12, clearly shows that the aircraft, when 
passing through a microburst, depending on the location of the flight 
path, will initially encounter increasing head winds, which then decrease 
and finally reverse direction to become a tail wind (see for example 
flight path #1 on Figure 2.7, page 13). Thus, the theoretical curves 
characterize the thunderstorm winds. 

To test the hypothesis that a varying wind having a frequency near 
that of the aircraft phugoid frequency could indeed cause the high 
amplitude and loss of control illustrated by some of the analytical 
models earlier, a number of computer runs and flight simulator tests 
were conducted. The frequencies of the sinusoids were 1, 1-1/2, and 2 
times the phugoid frequency of the aircraft type under study. 

4.4.3 Flight Path Deterioration Parameters 

In order to assess the potential severity of a wind shear hazard 
existing along a flight path, a quantitative parameter is needed to 
describe the response of the aircraft/pilot system. Variations of a 
parameter, referred to as a flight path deterioration parameter (FPDP) , 
was therefore defined. 

Table 4.2 gives expressions for the FPDP proposed for flight path 
and airspeed deviations. These are based on the following logic. 
Normalized altitude, HP/HG, was chosen as a flight path deviation 
parameter where HP is the height of the aircraft above the ground and HG 
is the height of the glide slope above the ground. The airspeed devia- 
tion parameter chosen is simply the deviation of airspeed from the 
reference value. However, in both flight path and airspeed parameters 
the positive deviations are examined separately from the negative devia- 
tions to avoid the cancellation errors. Also, the root mean square 
value of the velocity and height departure from the reference value was 
studied. 
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TABLE 4.2. Flight Deterioration Parameters Used in Comparing Computed 
Versus Manned Flight Simulator Control Performance in 
Idealized Thunderstorm Wind Shear. 


1 . AH 


L I 


(HP - HG)2dt 


where T|_ is the total landing time, HP is aircraft altitude, and 
HG is glide slope height. 


2a. 


GS+ = i 
n 


T 


HP 

HG 


dt 


where HP/HG above or on glide slope > 1 and T^ is the time above or 
on glide slope. T^/Tl is percentage~of time above or bn glide slope. 


Tm 


2b. GS- = 
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HP 

HG 


where HP/HG below glide slope < 1 and T^ is the time below glide 
slope. Tm/Ti_ is percentage of time below glide slope. 
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where is airspeed and V 
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(V - )dt 


is reference airspeed. 


for Va - Vao > 0 where T-j is the time airspeed is equal to or greater 
than reference airspeed. Ti/TL is percentage of time above or 
equal to reference airspeed. 



for Va - Vao ^ ^ where T|^ is the time airspeed is below reference 
airspeed. T|</Tl is percentage of time below reference airspeed. 
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4.4.4 Description of Test Plan 


In an effort to determine if the aircraft trajectory model simu- 
lates a real aircraft/pilot system, the results of the trajectory 
program were compared with a series of runs that were carried out in the 
B-727 simulator at the NASA Ames Research Center (NASA/Ames). The 
simulated aircraft were flown through the three different wave-form wind 
models characteristic of the thunderstorm downburst cell environment. 

The simulator runs were designed to test aircraft/pilot response to 
longitudinal and vertical wind waves of varying amplitudes and frequen- 
cies. The test plan used is given in Table 4.3. Twenty-seven computer 
runs and 19 manned flight simulator approaches were made. The computer 
and simulator test results are summarized below. Complete details of 
these runs are given in Turkel et al. [74]. 

Also, a series of manned flight simulations were conducted on a B- 
727 simulator at the United Airlines (UAL) Flight Training Center in 
Denver. All approaches were flown by a UAL simulator test pilot. 

Twelve B-727 ILS approaches were flown for a theoretical microburst 
single sine wave wind shear input. The head wind was first encountered 
at 430 m (1400 ft) AGL. The simulator phugoid frequency was 0.025 Hz or 
a period of 40 seconds. Wave amplitudes of 5, 10, 15, 20, 25, and 30 
m/s (10, 19, 29, 39, 49, and 58 kts) were flown. Eight of the 12 
approaches were flown at the 40-second period, while the remaining were 
flown at 10-, 20- , 80- , and 160-second periods each. 

Flight path trajectories measured in the NASA/Ames flight simulator 
are compared with values computed with the computer model in Figures 
4.10 through 4.14. Figure 4.10 compares computed and manned simulator 
flight path trajectories through a longitudinal sine wave of phugoid 
frequency and varying amplitude, i.e., 5.15, 10.3, and 15.45 m/s (10, 

20, and 30 kts). Figure 4.11 shows the same comparison for a longitu- 
dinal sine wave of 20.6 m/s (40 kts) amplitude and varying frequency, 
i.e., cop^, 1/2 ojp^, and 2 cop^. Similar comparison of computed versus 
manned flight simulator trajectories are given in Figures 4.12 and 4.13 
for a 1 - cosine downdraft wind shear. Figure 4.14 compares trajectories 
for a combination of longitudinal S-shaped and 1 - cosine downdraft winds. 
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TABLE 4.3. Test Plan for Simulator and Computer Runs 


Amplitude 


Aircraft trimmed for; 3.0“ glide slope 
70.0 m/s airspeed 
63,958 kg (140,000 lbs) 
gear down, flaps 30“ 

S-shape head wind to tail wind shear wave 


Full sine wave head wind to tail wind shear 


1 - cosine down draft 


Combinations: 

fS-shape, 0-5.15 m/s tail wind shear at 


jS-shape, 0-5.15 
\l - cosine, 5.1 


5 m/s (17 ft/s) downburst at w 


S-shape, 0-10.3 m/s tail wind shear at w 


|S-shape, 0-10.3 
ll - cosine, 10- 


3 m/s (34 ft/s) downburst at oj 


fS-shape, 0-20.6 m/s tail wind shear at a> 


:U.D HI/ 5 Lfl I I WlfiU aiicai au tUp^ 

ll - cosine, 20.6 m/s (68 ft/s) downburst at 

fS-shape, 0-20.6 m/s tail wind shear at 2 

20.6 m/s (68 ft/s) downburst at 2 



5.15 m/s 
10.30 m/s 
15.45 m/s 
15.45 m/s 
15.45 m/s 

5.15 m/s 
10.30 m/s 
15.45 m/s 
20.60 m/s 
20.60 m/s 
20.60 m/s 

5.15 m/s 
10.30 m/s 
15.45 m/s 
15.45 m/s 
15.45 m/s 


(10 kts) 
(20 kts) 
(30 kts) 


(17 ft/s) 
(34 ft/s) 
(51 ft/s) 
(51 ft/s) 
(51 ft/s) 




Computed Fliqht Simulator 


00 
98 
1.75 m/s 
0.70 m/s 


4.01 m/s 


3.57 m/s 


— Computed FI ight Simulator 


.50 

.73 

10.8 m/s 
7.1 m/s 


9.69 m/s 


4.20 m/s 


Sine Wave Amplitude * 5.15 m/s 


Frequency 


+ = 7.09 m/s 


4.78 m/s V 


3.19 m/s 


4.30 m/s 


9.0 m/s 


7.4 m/s 


6.55 m/s 


4.47 m/s 


Sine Wave Amplitude = 10.3 m/s 


2 


Frequency ■ 2 


Computed Flight Simulator 


1 .40 

GS+ - 

1.00 

0.79 

GS- - 

0.95 

8.99 m/s 

V+ - 

6.85 m/s 

-6.09 m/s 

V- - 

-5,'30 m/s 


Computed Flight Simulator 


GS+ = 1 .77 GS+ = 1 .02 

GS- = 0.94 GS- = 0.87 

V+ = 4.97 m/s V+ = 7.06 m/s 

V- = -3.78 m/s V- = -1.84 m/s 


Sine Wave Amplitude - 15.45 m/s 


Frequency * 1/2 


Figure 4.10 Comparison of computed 
and manned simulator 
flight path trajectories 
through a longitudinal 
sine wave of coph 
frequency. 


Figure 4.11 Comparison of computed 
and manned simulator 
flight path trajectories 
through a longitudinal 
sine wave of 20.6 m/s 
(40 kts) amplitude. 


























Figure 4.12 Comparison of computed 
and manned simulator 
flight path trajectories 
through a vertical 1 - 
cosine wave amplitude 
of 03 ^ frequency. 


Figure 4.13 Comparison of computed 
and manned simulator 
flight path trajectories 
through a vertical 1 - 
cosine wave of 15.45 m/s 
(30 kts) amplitude. 
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Computed Flight Simulator 




(b) 


— Computed FI ight Simulator 



Figure 4.14 Comparison of computed and manned simulator flight path 
trajectories through a combination longitudinal S-shaped 
and vertical 1 - cosine wave of 20.6 m/s (40 kts) 
amplitude. 
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An cases shown in Figures 4.10 through 4.14 were run for a modeled 
B-727 aircraft trimmed (flaps 30°) for a 3° flight path angle with an 
approach airspeed of 70 m/s (136 kts) and an angle of attack of 6.2°. 

All wind profiles are encountered at x = 0. The flight path deteriora- 
tion parameters pertaining to each case are given on the respective 
flight path trajectory plot. 

The computer model control system utilized thrust to control 
airspeed and elevator deflection to control flight path angle. It 
should be noted that fixed gains are used in the formulation of this 
model to represent an initial effort to model pilot control response to 
wind shear profiles and to simulate engine response characteristics. A 
discussion of the individual flight paths is given in Turkel et al. [74] 
and Frost et al . [72] . 

In terms of aircraft/pilot response, the computer model compares 
well with the simulator for the full sine waves, S-shaped waves, down- 
bursts, and combinations. However, some discrepancies exist with 
regard to the degree of flight path and airspeed control between the 
computer model and the test pilot. Although the control logic for the 
model pilot is similar to the control strategy of the test pilot, the 
test pilot flew consistently better than the model pilot. This is due 
to the fixed gain structure of the computer model pilot. A real pilot 
does not behave in the rigid manner of a fixed-gain model. In reality, 
a pilot acts in a variable gain decision-making process, which is 
probably not adequately included in the simplified models used in the 
study. This fixed-gain structure of the model allows for lower pilot 
damping of the flight path and airspeed oscillations induced by encoun- 
ters with wave disturbances. The test pilot is clearly of better skill 
than the computer model pilot. In addition, the test pilot had the 
opportunity to "learn" the types of profiles he was flying during the 
tests. 

4.4.5 Results of Flight Path Deterioration Parameters (FPDP) 

The airspeed deterioration parameters V+ and V-, calculated from 
the computer simulations, increase with increasing longitudinal wave 
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amplitude for waves at the phugoid frequency as shown in Figures 4.15 
and 4.16. Also shown is a comparison with the UAL flight simulator 
studies. The UAL flights were carried out only with sine waves and the 
comparison is thus limited. The largest V- values were attained for the 
2 ojpi^ waves, whereas the V+ value of the 20.6 m/s sine wave at 2 is 
not as large as that of the phugoid frequency wave. The values of V+ 
and V- from the computer simulations were smallest for the 1/2 ojp^ S- 
shaped and sine waves. In comparison with the computer results, similar 
trends are noted for the airspeed deterioration parameters, with the 
simulator runs for the S-shaped waves, and for V+ values of the sine 
waves. However, in the case of the sine waves, the simulator values of 
V- tend to be inconsistent. This is possibly due to the pilot "learning" 
the profiles and "fine tuning" his control procedures. The UAL manned 
simulator data is consistent in trend but considerably lower. The 
reason for this is believed to be due to the simulated low-frequency 
response in the phugoid range being overly damped in the training 
simulator. 


Control difficulty was encountered by the computer model and test 
pilot in flight through downbursts (particularly the 1/2 cop^ wave) and 
for the combination S-shaped longitudinal waves and 1 - cosine down- 
bursts. For the downbursts, the computer and the simulator airspeed 
deterioration parameters (V+ and V-) were lowest for the 1/2 ojp^ wave 
which, however, caused the largest deviation in flight path. This is 
reasonable since the long wave downburst does not have a pronounced 
effect on airspeed deviation but instead causes the aircraft to descend 
below the glide slope with the steadily descending air mass. The air- 
craft remains in this long wave for 76 seconds. Therefore, airspeed dete- 
rioration is probably not a meaningful warning parameter for application 
to downbursts. It may be noted that the glide slope deviation parameter 
GS- for downbursts shows very low values corresponding to large descent 
below the glide slope for the 1/2 cop^ waves. 

For the combination S-shaped longitudinal waves and 1 - cosine 
downbursts, the decreasing airspeed and the descending air mass forcing 
the aircraft below the glide slope presented the most difficulty for the 
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a) V- values versus sine wave amplitude and 
frequency. 



b) V+ values versus sine wave amplitude and 
frequency. 


Figure 4.15 Airspeed deviation parameters for longitudinal sine wave 
winds . 
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a) V+ values versus S-shape wave amplitude and 
frequency. 



b) V- values versus S-shape wave amplitude and 
frequency . 


Figure 4.16 Airspeed deviation parameters for longitudinal S-shaped 
wave winds. 
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simulated pilot and test pilot. Results show that the largest V- 
values for airspeed deviation correspond to the worst control cases and 
for the combined longitudinal and downdraft cases appears to provide a 
meaningful warning of hazardous wind conditions. 

Thus there is not a clear indication of which FPDP is most useful 
or whether a combination of parameters is required. The results of the 
computer and simulator runs through the longitudinal S-shaped waves and 
sine waves indicate a need for further studies to examine the effects of 
large longitudinal wind gradients due to large amplitude waves at short 
wavelengths. A parametric analysis of a broader range of wave frequen- 
cies must be carried out to determine the bandwidth which is most 
hazardous to aircraft operations. 

Control difficulties were noted with 1 - cosine downbursts, particu- 
larly the long duration wave at large amplitude and the strong down- 
bursts combined with longitudinal shear. Frost et al . [72] concluded 
that the combination longitudinal 5-shaped and vertical 1 - cosine wind 
profile is the most realistic profile of a downburst cell wind field in 
the vicinity of the ground. However, the high amplitude of the down- 
bursts studied may not be realistic close to the ground since the 
vertical wind component must approach zero there. New measurements of 
wind shear [32,39] will help provide meaningful magnitudes of the 
downburst. 

It is noted that because of the interrelationship of the longitudi- 
nal and vertical winds, indicated by the study, it is not clear that 
measurement of only the longitudinal wind component along the flight 
path for detection and warning of hazardous wind shear will be suffi- 
cient. This has severe ramifications since the vertical component is 
much more difficult to measure operationally. Proposed airborne and 
ground-based systems for detecting and warning of wind shear are 
discussed in the next section. All of these depend on measurements of 
only the longitudinal component. 
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5.0 DETECTION AND WARNING SYSTEMS 


5.1 Airborne Aids for Coping with Low-Level Wind Shear 

5.1.1 FAA Flight Tests for Airborne Aids 

Foy [19] reports a series of piloted flight simulation studies 
supported by analytical and experimental analyses of airplane response 
to wind shear and the meteorological phenomena producing low-level 
shear. Approach and landing tests were run under different conditions 
of visibility with different levels of approach instrumentation (full 
ILS and localizer only), and with wide-bodied and nonwide-bodied jet 
transports. The manned flight simulation experiments were run with a 
significantly large number of experienced pilots. 

A major conclusion over all the tests was that conventional (base- 
line) approach-management techniques, based on attempts to maintain a 
stabilized indicated airspeed from glide slope capture to the flare, are 
not effective in coping with the more severe (e.g., frontal and thunder- 
storm) wind shear encounters. Tests to develop improved approach 
management techniques considered both acceleration augmentation and the 
use of ground speed information. The results of these tests show that 
ground speed is particularly important. Although several potential 
solutions to the wind shear problem were indicated from the tests, the 
modified flight director with acceleration margin go-around indicator 
MFD/aA system performed well enough and ranked high enough in accept- 
ability to be recommended as a solution to the wind shear problem on 
approach and landing. 

The MFD/aA system contains the following combination of command 
i nformation: 

5. 1.1.1 Modified Flight D i rector . The modified flight direction 
includes improved flight-director control laws that incorporate accelera- 
tion augmentation to aid in coping with wind shear on approach and 
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landing [75], In comparison with the standard or baseline flight 
director commands, the modified steering control laws exhibit quickened 
response to changing wind and other transients. The modified flight 
director also has a modified speed command, driving the fast/slow "bug," 
that uses acceleration augmentation and wind shear compensation to 
improve speed control. For approach and landing, the pilot's speed 
control task is aided by supplying a speed-error indication on the 
fast/slow scale of the flight director. A basic assumption of the 
system is that a measurement of ground speed (GNS) is available in the 
airplane. 

5. 1.1. 2 Acceleration Margin . Acceleration margin, aA, is an 
analog quantity designed by FAA to indicate when the airplane is getting 
into a hazardous situation with respect to longitudinal wind shear. 


Acceleration margin is computed by: 


AA = - [-WD]fi/H 

(5.1) 

WD = (TAS - GNS) - 

(5.2) 

where 



A = acceleration capability of the airplane in level flight or in 
^ approach configuration (kts/s) 

WX . = wind component at the ground and along the runway (head wind 
^ is positive) (kts) 

TAS = true airspeed of the airplane (kts) 

GNS = ground speed of the airplane (kts) 

H = altitude of airplane center of gravity above ground; altitude 
is positive when measured upward (ft) 

A = rate of change of altitude with time; positive up (ft/s). 

In this case, A is a constant for the approach and will depend on the 
cap 

selected approach speed, the flap setting, the maximum engine thrust 
available, the drag, the aircraft weight, and the air density. For 
instance, values for the DC^-IO at 158,800 kg (350,000 lbs), 50° flaps, 
nominal approach speed, gear down, are: 
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2 

sea level, standard day, 0,86 m/s (1.7 kts/s) 

9000 ft, standard day, 0.51 m/s^ (1.0 kt/s) 

The term (TAS - GNS) is approximately the longitudinal wind velocity at 
the airplane (head wind positive), WD is thus the wind difference, or 
estimated wind shear, i.e., the difference in wind between the air- 
plane's present position and ground; a decreasing head wind is a positive 
difference. The magnitude of H/l^ is the expected time in seconds to 
reach the ground, and is negative for descent. Thus, the term 
[-WD]R/H is the expected acceleration demand due to longitudinal wind 
shear, with a decreasing head wind and a descending aircraft giving a 
positive demand. If the demand equals or exceeds A^^^p, aA becomes zero 
or negative and the situation is potentially hazardous. 

Tests with this system showed that the condition aA less than or 
equal to zero, if used as a criteria for advising a go-around, produced 
too many nuisance alarms. The algorithm was augmented with the differ- 
ence, DA, between the wind change and the airspeed pad given by: 

DA = WD - (IAS - V ) (5.3) 

^ app^ 

where 

IAS = indicated airspeed (kts) 

V = selected approach speed (kts) 

The go-around advisory is implemented according to Figure 5.1. The 
switches are closed when the indicated condition is true. The effect is 
to inhibit the go-around advisory if either the wind difference (decreas- 
ing head wind) is less than 12.9 m/s (25 kts) or the wind difference is 
no more than 4.1 m/s (8 kts) greater than the airspeed pad. The par- 
ticular values 4.1 and 12.9 m/s (8 and 25 kts) were chosen empirically. 

5. 1.1.3 Modifi ed Go-Around Guidance . The modified go-around 
guidance, intended to provide a pitch steering control law for use in 
wind shear, is based on the following rationale: 
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o 



WD>25 



Turn on 

Go-Around 

Light 


AA<0 


Figure 5.1 Go-around advisory augmentation algorithm [19]. 


0 The dominating requirement during go-around is terrain 
avoidance and obstacle clearance. After the initial 
pitch-up maneuver, it is assumed that flying a nominal 
positive flight path angle will result in a safe 
go-around. 

0 The pitch attitude required to maintain a flight path is 
dependent on the prevailing wind. The steering-control 
law should contain compensation for this effect. 

0 If there is severe wind shear or some other condition 
such that the aircraft cannot maintain the nominal flight 
path angle, the aircraft will be flown at or above a 
minimum airspeed at a commensurate maximum pitch attitude. 

Vertical speed, H, and ground speed, GNS, inputs were used to 
compute flight path angle, y. Flight path angle and angle of attack, a, 
were input into the computation of the pitch steering signal, A. This 
signal and the pitch rate term, e, are the controlling terms for 
damping as long as the airspeed remains high. When airspeed drops to or 
below the stall value, a minimum function selector chooses the IAS - 
Vstaii input, which results in a pitch-down command to gain airspeed. 

The reference flight path angle, y^^, and angle of attack, were 

chosen empirically. 

With the modified go-around method, the pilot advances the throttles 
to give full thrust immediately after deciding to go around. He is then 
not using the F/S indicator on the flight director for the thrust 
control. Therefore, to provide additional information, the F/S signal 
was modified so that the F/S displayed an approximation to angle of 
attack error. 
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The MFD/aA system, which showed a significant performance improve- 
ment over baseline in the wind shear studies [19], requires instrumenta- 
tion to measure certain aircraft variables and wind components that are 
not available in many current aircraft. Of the quantities that are 
usually not available or are not measured adequately, the most important 
is ground speed, altitude above the runway, and rate of change of alti- 
tude. Additionally, there is a firm requirement for accurate knowledge 
of the winds on the runway; the along-runway component is needed by 
algorithms such as the acceleration margin and the crosswind component 
to enable the pilot to anticipate his lateral control action. 

Foy's tests showed importantly that there are realistic wind shear 
conditions that can occur on takeoff which exceed the aerodynamic lift 
and thrust capability of the airplane. An attempt to make a normal 
takeoff in such a situation, even when aided by a minimum height loss 
pitch-steering algorithm, cannot be handled by pilot action. The most 
appropriate recourse found in the study is: (1) not to attempt to 

takeoff at all, (2) to take off in a different direction, or (3) to 
prolong the takeoff roll so that rotation will lift the airplane off 
with 10.3 m/s (20 kts) or more of excess airspeed. Any of these actions, 
in practice, requires advance notice (that is, prior to starting the 
takeoff roll) of the wind shear condition. 

5.1.2 Safe Flight Instrument 

A self-monitoring wind shear warning system has also been developed 
by Safe Flight Instrument Corporation (Stein [76]; Greene [77]). This 
system is designed to sense and integrate horizontal and vertical, or 
downdraft, wind shear components providing the pilot of an aircraft on 
approach a timely warning to initiate a go-around. 

The wind shear monitoring system computes the thrust required to 
maintain the desired glide path when a downdraft is encountered on 
approach. The thrust required in g's is equivalent to the angular 
displacement from that glide path when the actual (or potential) devia- 
tion is measured in radians. This displacement, termed downdraft drift 
angle (DDA), is a function of the ratio of the velocity of the descending 
air to the aircraft's speed; 
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(5.4) 


DDA = W^/V^ 
where 

DDA = downdraft drift angle 
= vertical wind 
Vg = airspeed 

The effect on the airplane's landing profile due to a change in the 
head wind (tail wind) component due to wind shear may be described by 
the rate of change of ground speed (inertial acceleration) required to 
maintain a constant lift condition (airspeed acceleration = 0). The 
magnitude of a horizontal wind shear is predicted by Greene [77] with 
the following formula: 

WSx = (V^ - x)/(H/H) (5.5) 

where 

WS = horizontal wind shear 

A 

Vg = airspeed 
X = ground speed 

H = height or bandwidth of the shear layer 
^ = vertical velocity 

Figure 5.2 from Greene [77] shows the functional block diagram of 
Safe Flight Instrument Corporation's wind shear computer. The computer 
resolves the two orthogonal vectors of a wind shear encounter and 
provides meter output and threshold alert indication of that encounter. 
The two vectors are called DDA and horizontal wind shear (HWS) . 

HWS is derived by subtracting longitudinal acceleration from 
airspeed rate. Airspeed rate is obtained by taking an airspeed analog 
from the airspeed indicator, or air data computer, and passing it through 
a high-pass filter. Longitudinal acceleration is sensed by a computer 
integral accelerometer, the output of which has been summed with a pitch 
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Figure 5.2 System block diagram of Safe Flight Instrument Corporation’s 
wind shear computer [76] . 
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attitude reference from a vertical gyro signal to correct for the accel- 
eration component due to pitch (g sin e). The summed acceleration and 
pitch signals are fed through a low-pass filter, the output of which is 
summed with the airspeed rate signal to comprise horizontal wind shear. 

The vertical component of DDA is developed through the comparison 
of measured normal acceleration with calculated glide path maneuvering 
load. Flight path angle is determined by subtracting the pitch attitude 
signal from the angle of attack analog as sensed by the stall warning 
airflow sensor. This is then introduced into a high-pass filter and 
then a multiplier to which the airspeed signal has been applied. Thus 
the flight path angle rate, corrected for airspeed, provides the computed 
maneuvering load term. This term is compared in a summing junction to 
the output of a computer integral normal accelerometer. A failure to 
match is the indication of an acceleration due to downdraft. If this is 
the case, this acceleration when integrated, is the vertical wind 
velocity and is further divided by the airspeed signal to compute DDA. 

The DDA and horizontal wind shear signals are combined and the 
summed output passes through a low-pass filter forming the output 
signal. This signal is fed to a comparator which provides a latched 
ground output signal (for a warning device) and a meter output. The 
warning output is set at a threshold of -0.67 m/s (-3 kts/s) hori- 
zontal wind shear DDA of -0.15 rad or any combination which would total 
an equivalent signal level. 

5.1.3 Bliss's Aircraft Control System for Hind Shear 

Bliss [53] questions whether acceleration augmentation and quicken- 
ing of pitch steer commands are sufficient to solve the wind shear 
problem. He believes that modified flight directors, as used in the FAA 
wind shear experiments, using the conventional IAS parameters for the 
approach speed, will result in the same hazardous ground speed values 
close to the ground and will produce the same results as exist today. 
According to Bliss, a flight director utilizing a totally computerized 
inertial vector, vertically, laterally, and longitudinally (speed 
vector), wherein the ground speed is integrated properly with the 
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indicated airspeed, is needed to resolve the wind shear problem. 

Bliss states that the minimum standard instrumentation for certification 
is as follows: 

An analog groun d speed instrument mounted in close 
proximity to, or combined on the same instrument with, 
an analog airspeed instrument. 

2. A system of three indexes: 

An airspeed target index selectively adjustable 
by the pilot to indicate the normal minimum 
approach indicated airspeed value. 

b) A first ground speed index automatically pro- 
grammed to a ground speed value equal to the 
true airspeed value of whatever the IAS index 
is set on. (This then, becomes a zero wind 
ground speed index value ) . 

c) A second ground speed target index programmed 
to a ground speed value relative to the zero 
wind index, taking into account the surface 
head-wind/ tail -wind component on the runway. 

This index is the ground speed expected 
approaching the threshold, and it then becomes 
the minimum ground speed value for that approach. 

3. The use of two minimum speed values requires that they 
be automatically integrated to a third instrument which 
then becomes the primary speed instrument . The use of 
which eliminates the use of speed values slower than 
either the normal approach minimum airspeed or the 
normal approach minimum ground speed . (This can be a 
fast/slow instrumehTT] 

4. A tail wind warning system variably programmed with alti- 
tude, which calls the pilot's attention to the excess 
ground speed existing on the approach when it is not 
possible for the aircraft to decelerate inertially to 
normal values before reaching the landing point. 

5. An excess head wind warning (programmed much the same as 
the tail wind warning) of values of excess IAS variably 
programmed with altitude, to warn the pilot when his 
elevator control authority will be limited after the loss 
of airspeed results in normal airspeed. This warning may 
contain a limiter when the airplane is trimmed to a 
nose-up trim with the excess airspeed so that after the 
airspeed loss, the airplane will be in an acceptable 
trim condition. 
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6. Airports served by Part 121 carriers must by required to 
give surface wind information in the landing area for 
landing traffic , and in the vicinity of the departure 
end of the runway for aircraft taking off . They must 
also have remote wind sensors located at the highest 
elevation possible on any obstruction requiring unusual 
climb-out procedures'. 

7. Aircraft operating under Part 121 which are equipped 
with INS must have a recording on the flight data 
recorder of a ground speed parameter. 

8. All Part 121 aircraft must have an on-board ground speed 
detection system capable of an accuracy of less than 2 
percent error and a ground speed tracking error of less 
than 1 . 5 sec. 

9. All certification of auto-land systems should be canceled 
until they are modified to the standards provided by this 
airspeed/ground speed system, including the full pilot 
monitoring instrumentation. 

10. For the proper solution to the wind shear problem in all 
aircraft, a standard means for providing ground speed 
must be adopted . The least expensive (even light 
trainers may use it), may be an airborne Doppler-type 
system with a compatible ground-based transponder. The 
ground-based transponder can be located at the inter- 
section of two or more runways for use in any appropri- 
ate direction. 

5.1.4 Advantages and Disadvantages of Airborne Systems 

The preceding discussion relates to airborne systems. There are 
several important advantages to an airborne system: 

1. Each aircraft properly equipped with an airborne system 
carries its protection wherever it flies. Thus, a 
ground-based system is not required at each airport. 

2. The system allows the pilot to monitor the changing 
longitudinal wind shear conditions in a quantitative 
manner, that combines shear and aircraft performance in 
a meaningful way. 

3. Some advanced indication, even if only a few seconds, is 
given to the pilot, so speed banking and/or a go-around 
can be attempted. 

The disadvantages, however, include: 
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7. The system requires a sophisticated ground speed measuring 
system for the aircraft. In the U.S. Civil Airline Fleet, 
essentially only the new wide-bodies transporters currently 
have such a capability. The larger number of smaller jets 
do not have such ground speed measurement capabilities. 

It is unclear that an acceptably accurate, inexpensive 
system can be developed for these aircraft, since inertial 
or other high-resolution navigation systems are likely 
required for the measurement. Finally, the requirement 
for ground speed measurements would be even more difficult 
to achieve for those aircraft in the general aviation 
fleet susceptible to wind shear (i.e., light-business jet 
transports) . 

2. The acceleration margin system requires an airborne wind 
measurement system, clearly requiring an inertial -type 
measurement for sufficient accuracy and resolution. 

3. The system requires notice of the runway threshold wind 

to be. given to the pilot. For microburst events and other 
small wind shears which can occur very rapidly, a few 
seconds delay in updating runway wind can seriously 
hamper the system's effectiveness. A telemetering of 
runway wind probably is needed which thus results in the 
requirement of equipment at each airport, and hence 
removing part of advantage #1 listed earlier. 

4. The system makes the assumption that the longitudinal 
wind shear component is sufficient to determine the 
threat. As discussed earlier, some uncertainty remains 
concerning this point. 

5. Perhaps the most serious limitation lies in the fact that 
during takeoff, wind shear so severe that a suitable 
acceleration margin is unavailable for aircraft survival, 
can be readily encountered. Also during approach an 
aircraft must enter a dangerous wind shear condition 
before having the data to make corrective action. 

6. In the flight simulator testing of this system, more 
realistic wind shear profiles need to be tested. 

7. Using the acceleration margin technique, the presence of 
a phugoidal instability forcing in the wind shear is not 
considered in A^^p- Thus, further, theoretical and flight 
simulator testing of the concept is required. 

Despite the shortcomings inherent in the airborne system, it 
probably provides the best detection/warning capability to date for an 
aircraft in flight, and undoubtedly aids the pilot by providing up-front 
data to aid in traversing severe wind shear conditions. Although the 
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system may not be effective in some situations, for many others it may 
clearly save the aircraft. 

5.2 Ground-Based Wind Shear Detection 
and Warning Systems 

In the past few years, a number of ground-based wind shear 
detection/warning systems have been proposed and some tested. Notably 
among these is the low-level wind shear alert system (LLWSAS), the 
thunderstorm gust front detection systems based on combinations of wind 
and pressure sensors, the acoustic Doppler system, the laser system, 
and the pulsed microwave Doppler radar system. 

5.2.1 Low-Level Wind Shear Alert Systems (LLWSAS) 

The LLWSAS is an operational FAA near-term solution to the wind 
shear hazard. The LLWSAS detects the presence of wind shear in the 
vicinity of the airport at the surface. Plans to install 51 more of 
these units at major airports within the United States are underway. To 
date, 58 systems have been installed. 

The system consists of an airport-centered array of six anemometers 
clustered at approximately 3-km (2 mi) spacing with a reference sensor 
located near the geographic center of the airport. The data are tele- 
metered to a master station in the control tower and processed by a 
minicomputer. If the LLWSAS computer senses a vector difference of 15 
kts or more between the mid-field and perimeter winds, it activates an 
aural alarm and a display screen in the control tower. A warning is 
then transmitted to the pilot by an air traffic controller. 

The LLWSAS system, however, cannot guarantee protection in all 
cases’. On August 22, 1979, an Eastern Airlines B-727 on approach to 
William B. Hartsfiel d/Atlanta International Airport dropped suddenly 
from 750 to 375 ft above ground level in a strong shear despite the 
flight crew's immediate decision to execute a missed approach. The 
LLWSAS on the airport remained mute. The fact that the system does 
not measure the wind shear at a height above the surface, where the 
actual aircraft problem exists, is not just a limitation; it creates 
the potential for false security, which does not exist. 
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The system, moreover, was designed to detect large horizontal wind 
shears that move across the airport, as seen in surface wind data. 

Thus, the system is suited for cold frontal passage and thunderstorm 
gust fronts but is not well suited to detect smaller scale phenonema 
such as the outflow portion of a microburst. It is equally apparent 
that the LLWSAS is unable to detect the downdraft associated with 
microbursts or other forms of vertical winds. 

5.2.2 Pressure Jump System 

This system is based upon the characteristic pressure jump that 
proceeds frontal wind shear. The system comprises a large array of 
pressure jump detectors distributed in a dense pattern around the 
ai rport . 

Although the system has proven to be rather successful in gust 
front detection, false alarms resulting from turbulent wind gusts and 
certain technical difficulties have caused delays in implementation of 
the system. 

5.2.3 Acoustic Doppler System 

The acoustic Doppler system determine wind speed and direction by 
measuring frequency shift (Doppler effects) in signals reflected by the 
atmosphere. The system was found to be expensive and unable to operate 
under heavy precipitation and in zones of noise created by aircraft. 

5.2.4 Laser Systems 

The laser system scans directly over the sensor using a continuous 
wave laser. This system does not have the range required to scan the 
glide slope and takeoff flight path to detect wind shear. There is a 
possibility that this capability may be available in the future using 
a pulse Doppler laser technique. 

5.2.5 Pulse Microwave Doppler Radar 

McCarthy, et al. [78,79], Wilson et al . [80,81], Fujita and Wakimoto 
[25], Offi et al. [82], and Strauch [83] have all demonstrated the utility 
of ground-based pulsed microwave Doppler radar to measure low-level wind 
shear events. McCarthy et al . [78] used a NSSL Doppler radar to measure 
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wind data along the precision approach path to the Norman, Oklahoma, 
airport and verified measurements with two instrumented NCAR aircraft. 
Similar results are reported by Off i , et al . [82]. Comparison of 
Doppler- radar-measured winds with that measured by an aircraft are shown 
in Figure 5.3. Frost and McCarthy have proposed a detection and warning 
system which utilizes ground-based Doppler-measured wind data to predict 
aircraft performance. 

The proposed operational detection and warning system operates on 
the following principles: The wind speed profile is measured in real 

time with a Doppler radar looking along the flight path. The Doppler 
radar takes a wind measurement in 150-m (500 ft) steps (approximately 
every 2 seconds of an aircraft trajectory at 72 m/s (140 kts) approach 
speed). The wind data can be transmitted to either the approaching 
aircraft or to the air traffic controller. However, more optimum is a 
minicomputer or microcomputer slaved to the Doppler which applies air- 
craft response functions to the wind profiles for specific aircraft type 
and simulates aircraft trajectories. The flight path deterioration 
parameter based on the techniques described earlier (see Section 4.0) is 
determined in aeal time. An excessive value of the parameter triggers a 
warning alert. Figure 5.4 conceptually illustrates the technique. Some 
questions which remain to be resolved prior to developing an operational 
system are: (1) Is the longitudinal wind speed component more signifi- 

cant than the vertical component? (2) What is the definition of a 
meaningful flight path deterioration parameter? (3) What is the most 
complete and computationally efficient flight trajectory computer 
program for real-time application to computing flight deterioration 
parameters? 

The advantages of this concept are: (1) It quantifies the wind 

shear in terms of actual aircraft performance; (2) it provides a warning 
to an aircraft prior to the aircraft beginning the approach, as needed 
with the airborne systems; (3) the Doppler directly measures the wind 
along the glide slope and is not limited to the surface measurements; 

(4) it provides a numerical classification as to aircraft type (flight 
path deterioration parameter); (5) provides service for all sections of 
aviation, i.e., general aviation, corporate aviation, as well as 
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Figure 5.3 Comparison of aircraft (solid line) longitudinal wind 

and Lagrangian Doppler velocity (dashed line) as a function 
of time, for May 16, 1979, as part of SESAME 1979 [79]. 
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commercial airliners; (6) the system provides capabilities for both 
ground-based and airborne displays (data uplink); and (7) the system is 
an all-weather system. 

Two advantages of this system were strongly supported by the air 
traffic control committee at the Third Annual Workshop on Meteorological 
and Environmental Inputs to Aviation Systems [84]. These are: (1) A 

ground-based detection system must be able to detect wind shear along 
the approach to and departure from the runway and at an altitude to 
support the en route air traffic control system; and (2) wind shear 
intensity should be reduced to a numerical value which the pilot can use 
to determine if the intensity of the system is too great for his type of 
aircraft to penetrate (which currently is operationally undefined). 

Some disadvantages of the sytem include: (1) The system best 

measures the radial or longitudinal component along the intended approach 
path, the vertical component or downdraft cannot be measured directly in 
the current system; and (2) to utilize this system, each airport must be 
equipped with a Doppler radar, which can be a substantial expense. 

5.3 Current Status of Low-Level Wind Shear 
Detection and Warning Systems 

Although all the reported wind shear detection and warning systems 
have merit, no one system has proven to be fully adequate for fail-safe 
detection of low-level wind shear. Many of the systems are preliminary 
solutions which have been partially implemented without a thorough 
understanding of the nature of the problem. 

As noted, the LLWSAS and pressure jump systems do not measure the 
environment above the surface in which the aircraft may encounter wind 
shear. Moreover, they probably do not provide protection for the small- 
scale microburst-type wind shears. A relatively negative consequence of 
these two ground systems may be that they provide confidence for the 
pilot and controller in a system that may be less than adequate for 
certain dangerous situations. Moreover, the designers of the system may 
well understand the limitations but the users may not. 
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The various airborne systems are extremely useful but are not fail- 
safe. Based upon the method of storing kinetic energy to overcome 
sudden airspeed losses occurs only in the case when enough energy can be 
banked to accelerate the aircraft faster than the wind is decaying. 
Obviously, this does not work for takeoff. Moreover, the detailed 
flight simulation studies of this system may have some inherent disad- 
vantages. The wind shear models utilized in perfecting the MDF/aA, 
i.e., modified flight director acceleration margin system, are incomplete 
wind shear models. These models contained neither the lateral variations 
in wind, the appropriate turbulence intensity and distribution over the 
aircraft, nor include the very localized intense short-duration micro- 
burst which has been clearly identified from radar Doppler measurements. 
Moreover, the question remains as to how significant is the.phugoid 
oscillation of the aircraft. Aerodynamicists and pilots frequently 
point out that the phugoid oscillation is of such low frequency that it 
can easily be controlled. They have not considered the fact that 
forcing the aircraft with a forcing function, having the frequency of 
the phugoid, can appreciably augment the difficulty to control the 
subsequent motion. In turn, many flight simulators do not appropriately 
model the phugoid oscillations. If these oscillations are not appro- 
priately modeled by the simulator, then the performance of the aircraft 
will be quite different in a wind shear forcing the aircraft at this 
frequency than the flight simulator would demonstrate. 

Another question associated with the airborne system is: Can the 

airlines absorb the high cost of implementing the ground-speed and 
aircraft-speed measuring systems? Moreover, there must be some method 
for providing in real time the runway threshold winds to the aircraft. 
Finally, the very important question which must be resolved is whether 
a warning and detection system, either ground based or airborne, is 
adequate if only a measurement of the longitudinal wind speed is uti- 
lized, i.e., how significant is the vertical wind speed component in 
creating hazardous flight conditions? 

The ground-based pulse microwave Doppler can provide a great deal 
of help to the wind shear detection and warning requirements. It 
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provides the high-resolution detail of low-level shear, which can be 
processed to give predicted aircraft response, for either potential 
approach to landing and takeoff modes, without the aircraft actually 
entering the expected hazardous airspace. Its capability, however, has 
not been fully tested. Studies need to be carried out to determine 
whether it would ultimately be cost effective to the public to have a 
dedicated Doppler for the airport environment. 

There is also the problem of the basic theoretical concepts of 
aircraft performance in wind shear conditions. Practically all aircraft 
analyses are based on steady or zero wind conditions. A better under- 
standing of the ability of various aircraft to survive wind shear is 
necessary. There is a strong probability that several general aviation 
aircraft accidents, where flight data records are not available, have 
occurred due to wind shear and have gone undetected. Wind shear models 
typically utilized are two-dimensional steady-state models. Most all 
aircraft analyses have utilized three-degrees-of-freedom systems of 
equations. The microburst is clearly not a simple two-dimensional model 
but highly three-dimensional and time-dependent as well. Wind shear 
data input to numerical simulation models and to flight simulators must 
therefore be improved. Improvement in such models is particularly 
important with the airlines moving toward nearly 100 percent reliance on 
flight simulators for their training and proficiency needs. The flight 
procedures for pilots when encountering wind shear must be totally 
developed in manned flight simulators. Severe wind shear would be 
encountered by most pilots once, if at all, in a lifetime. However, if 
a wind shear encounter under realistic conditions is mandatory during 
flight simulator training, the pilot will have a much better concept of 
wind shear and be less likely to take lightly any wind shear alert 
warnings he may receive. 
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6.0 CONCLUSIONS 


Based on the review of wind shear hazard studies, the following 
conclusions have been reached. 

1. Current mathematical wind shear models, fast becoming 
standards, are not three-dimensional and are based on a 
few highly smoothed data. Turbulence superimposed on 
the wind shear is artificial and does not simulate the 
extreme turbulence reported by Bliss during his approach 
to Kennedy prior to the Eastern 66 accident. 

Realistic time-dependent three-dimensional wind shear 
models based on complete data sets are needed to fully 
verify airborne warning and detection systems, to develop 
flight procedures, and to train flight crews. The NASA 
Gust Gradient and NCAR JAWS programs have provided some 
data sets, but they remain to be analyzed. 

2. Order of magnitude analyses of the aircraft equations of 
motion show that horizontal wind shear terms generally 
produce the largest forces disrupting flight. These 
analyses suggest that values of horizontal wind shear 
smaller than 8 kts/100 ft, given in AC-20-57A for certi- 
fication of automatic control systems, can be critical. 
Values of wind shear should be specified as applied along 
the line of flight. Currently, vertical variation of 
wind speed is implied in AC-20-57A. 

3. Disagreement exists relative to the optimum flight proce- 
dures to employ when caught in wind shear. The argument 
of trading velocity down to stick-shaker speeds to enhance 
climb is opposed by the ALFA Airworthiness and Performance 
Committee who argues best climb performance occurs at 
minimum drag speed. The committee's premise is that 
flying at this speed leaves some excess kinetic energy or 
velocity to flair the aircraft at the last moment if 
impact is unavoidable. These arguments are primarily 
based on performance analyses using charts which are valid 
for 1-g flight conditions at constant indicated airspeeds. 
Dynamic analyses with realistic wind shear models are 
required to clearly resolve optimum flight procedure in 
severe wind shear. 

4. Simple mathematical studies of aircraft motion without 
control laws using linear as contrasted to logarithmic 
vertical wind speed profiles show conflicting results. 
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The reason for this is that the initial or trimmed flight 
conditions remain constant without control input. For 
the linear wind speed profile, the wind shear term is 
also constant and the aircraft remains in trim relative 
to the magnitude of the wind shear. For the logarithmic 
profile the wind shear term changes continuously along 
the flight path and the aircraft is thus always out of 
trim relative to the wind shear. 

5. Linear stability analysis clearly indicates that wind 
shear strongly affects the phugoid stability of the air- 
craft. This is further verified by nonlinear analysis 
which shows strong amplification of the phugoid oscilla- 
tion in typical thunderstorm-type wind shear. Pilot 
models and automatic control laws can, in general, cope 
with these oscillations although they may become uncon- 
trollable if the simulated pilot's skills are low or if 
the control laws lack sophistication. 

6. Computer and manned flight simulator studies of aircraft 
performance in sine and half-sine wave longitudinal winds 
and 1 - cosine downdrafts were carried out. The sinu- 
soidal wind profiles were applied along the flight path 
and represented a hypothetical representation of thunder- 
storm wind shear. The results of this study show that 
both the computer models and the manned flight simulators 
have most difficulty coping with combined longitudinal 
and downdraft wind shear. The second most difficult wind 
condition was the 1 - cosine downdraft used alone. These 
results suggest that a wind shear warning and detection 
system must measure the vertical wind component as well 
as the longitudinal component. Serious implications are 
inherent in this observation because the vertical wind 
component is much more difficult to measure than the 
longitudinal component. The magnitude of the downdraft 
velocities for the hypothetical wind shear models were 
chosen somewhat arbitrarily, however, and further study 
is required using realistic downdraft wind shear models 
to fully verify this conclusion. 

7. Six flight path deterioration parameters, FPDP, defined 
as a measure of the severity of a given wind shear condi- 
tion on aircraft performance were tested. Both computer 
analyses and manned flight simulator studies were carried 
out which showed general correlation between the magnitude 
of the FPDP and the quality of the computed and measured 
flight paths. In general, the FPDP defined as the root 
mean square difference between actual airspeed and refer- 
ence approach airspeed showed the best correlation with 
hazardous conditions for longitudinal wind shear. In 
turn, the FPDP defined as the difference of the actual 
flight path height minus the intended glide slope height 
divided by the aircraft's absolute altitude served as a 
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better measure of flight deterioration in downdrafts. 

All studies of the FPDP's were carried out using hypo- 
thetical sinusoidal wind shear models of different fre- 
quencies and amplitudes. Studies using realistic wind 
shear models based on measured data are needed to fully 
determine a realistic FPDP or to define an alternate 
measure of the severity of the wind shear to the aircraft 
performance. A quantitative value of the FPDP which can 
be computed in real time with a microcomputer "slaved" to 
a Doppler radar measuring the wind speed along the flight 
path is believed to promise the most effective operational 
system for warning of wind shear hazards. When the crit- 
ical value of the FPDP parameter is exceeded, a warning 
alarm would sound in the control tower and in the TRACON 
as well, if needed. In general, flight controllers prefer 
a numerical value of a warning parameter which can be 
used in the above fashion. 

8. Airborne systems developed to date measure only the longi- 
tudinal wind speed component and, in general, incorporate 
the concept of conserving energy for the situation when 
the wind shears from a head wind to a tail wind. These 
systems will not work during takeoff where maximum or 
essentially maximum thrust is already employed. Although 
the system has the advantage of being carried with the 
aircraft such that the warning and detection system is 
available regardless of where the approach or takeoff is 
made, it has the disadvantage that one must enter the 
hazardous airspace prior to the system providing any use- 
ful information. Additionally, the system requires a 
highly accurate ground speed measurement which is not 
generally available on the majority of commercial air 
carriers. In turn, most of the airborne systems have 
been developed and verified in manned flight simulators 
using incomplete wind shear profiles. There is clear 
evidence that wind shears can be encountered which are 

so severe that a suitable acceleration margin is unavail- 
able for the aircraft to survive the wind shear encounter. 

9. The current ground-based low-level wind shear alert 
system, LLWSAS, is only a near-term solution to the wind 
shear hazard. The LLWSAS system does not measure the 
environment above the surface in which the aircraft may 
encounter wind shear, and moreover increasing evidence 
illustrates the many severe wind shears are of sufficiently 
small scale, i.e., microburst-type wind shear, that they 
can occur directly over the airport and ao undetected by 
the LLWSAS. 

10. The ground-based pulse microwave Doppler promises to pro- 
vide the most effective wind shear detection and warning 
capability. The Doppler has been demonstrated to provide 
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high-resolution details of low-level wind shear which can 
be processed with a microcomputer to give predicted air- 
craft performance. Both the approach and takeoff modes 
can be handled without the aircraft actually entering a 
hazardous wind shear condition. The Doppler radar directly 
measures the wind along the glide slope and is not limited 
to a surface measurement such as the LLWSAS. Obviously, 
one Doppler radar can only measure the radial or longi- 
tudinal component of the wind along the intended approach 
path. If the vertical or downdraft wind speed component 
must be measured directly, two radars are required. The 
implementation of all major airports with one Doppler 
radar will be expensive but feasible; two radars may be 
prohibitive. Additional studies to fully develop the 
flight path deterioration parameter concept and to illus- 
trate that one Doppler radar per airport is sufficient 
to serve as a warning and detection system, are required. 

The hazard of wind shear to aviation operations is far from solved. 
The LLWSAS system currently installed at 58 airports may give too many 
false alarms and consequently causes complacency relative to wind shear 
situation. Airborne systems have not been implemented to any major 
extent and in turn may provide a sense of capability to cope with wind 
shear which is not real. 

The NCAR JAWS and NASA Gust Gradient field programs have provided 
the necessary data to make a quantum step forward in solving the wind 
shear problem. These data must be thoroughly analyzed, however, and 
appropriately formulated to allow development of effective warning and 
detection systems to provide mathematical wind shear models for flight 
crew training and to establish aircraft design criteria. Analysis of 
these data should proceed as rapidly as possible before further catas- 
trophes occur due to insidious wind shear lurking in the approach and 
takeoff air corridors of our major airport terminals. 
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APPENDICES 


APPENDIX A 


GENERAL EQUATIONS OF UNSTEADY MOTION 


Etkin [4] gives a complete development of the general equations 
of unsteady motion. However, the variation of wind velocity is not 
generally incorporated into the equations, i.e., a zero or constant wind 
is assumed. 

In this study, incorporation of the wind vector components into the 
governing equations is discussed. The set of equations is based on the 
assumption that the earth is a stationary plane in inertial space. This 
assumption is well justified for takeoff and landing problems. A coor- 
dinate system fixed at the earth thus becomes the inertial frame of 
reference, designated F^. The vehicle is assumed to be a rigid body 
having a plane of symmetry. 

In establishing the appropriate reference frame for computing the 
motion of the aircraft subject to a ground wind, we are particularly 
interested in an atmosphere-fixed reference frame, F^, since the aero- 
dynamic forces depend on the velocity of the vehicle relative to the 
local atmosphere. If the atmosphere is in uniform motion with velocity 
relative to the earth, then F^ moves relative to F^ with that velocity. 

Two other reference frames of interest are the air-trajectory refer- 
ence frame, F^^ (also called the wind-axis reference frame; this "wind" 
should not be confused with the atmospheric motion), and the body-fixed 
reference frame, Fg, or body-axis reference frame. The wind-axis refer- 
ence frame, F^, has the origin fixed to the vehicle, usually at the mass 
center, and the axis is directed along the velocity vector of the vehicle 
relative to the atmosphere, Thus, 

V = Vg - W 

where is the inertial velocity or the velocity of the vehicle relative 
to the fixed earth. The axis 0 z lies in the plane of symmetry of the 

W Vf 

vehicle. The frame F^ has angular velocity relative to the inertial 
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frame, F^, the components of which are conventionally designated by p^, 
and 

The body axes are in a body-fixed reference frame in a rigid body. 
Bodies with articulated control surfaces and/or elastic motions for which 
the body cannot be taken as rigid are not considered in this equation 
development. 

The origin of the body axes is usually the mass center of gravity, 

C. The plane of symmetry is generally taken as Cxz, with z directed down- 
ward. By convention, the components of angular velocity of the body-axis 
frame of reference, Fg, relative to F^ are designated p, q, and r and the 
components along the body axis of aircraft velocity relative to the atmo- 
sphere frame of reference, F^, are denoted by u, v, and w. 

On these assumptions, the classical six-degrees-of-freedom equations 
of motion (Equations 5, 8, 1, and 7 of Etkin [4] with the atmospheric 
wind effects included) become: 


Force Equations in Wind Axes, F|^: 

■^xw - ° - ""S 5''" ®w = 

- C + mg cos 6^ sin = mlSy^ + m[r^{V + 

- L + "ig cos 0^ cos + in[p„Wy„ - %(V + W^^)] 

Force Equations in Body Axes, Fgi 


(A. 2a) 
(A. 2b) 
(A. 2c) 


X - mg sin 0 = m(u + 1^1^) + m[q(w + W^) - r(v + W^)] (A. 3a) 

Y + mg cos 0 sin <J) = m(v + 1^^) + m[r(u + W^) - p(w + W^)] (A. 3b) 

Z + mg cos 0 cos 4> = m(w + 1^^) + m[p(v + W^) - q(u + W^)] (A. 3c) 

The components of the wind velocity vector are most frequently given in 

the earth frame of reference. The relationships between the earth compo- 
nents and those in the aircraft wind frame of reference, F^, are given by: 

“xw = “xE ®w ^ V ®w - ”zE ®w 
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% " '*’w ®w ""w ■ '*'w^ 

+ Wyg(sin 4.^ sin e^ sin + cos 4,^ cos ipj 

+ sin <f)^ cos 0 ^ (A. 4b) 

”zw = '^xe'"°® *w ®w *w "■'■'’ *w> 

+ Wy^(cos ((.^ sin e^ sin lO^ - sin 4 .^^ cos <(;^) 

+ cos <p^ cos 0^ (A. 4c) 


The velocity vector components in the body frame of reference are the 
same with the Euler angles (c{)^, 0^, 4»^) replaced by (<j), e, i{j). 

The total derivatives of the wind vector components are: 


'9W 1 


faw 1 

V 


f9W 1 

Y 


^ . t ^ 

aw 

Y 

[sx J 

Xp + 

E ^ 

A 

[9y J 


A 

[9Z J 


A 

[at J 

fsw 1 


f9W 1 


faw 1 


faw 1 

y 

[9X J 

Xp + 

'e 

y 

[9y J 


y 

az 


y 

[at J 


faWzl 


faw ' 


^aw ' 


'aw ' 

[3X J 

l^E^ 

[sy J 

^E 

F ^ 

lazi 

^E ^ 
F 

[atl 


(A, 5a) 


(A. 5b) 


(A. 5c) 


Additional equations are: 

Moment Equations in Body Axes : 

L = IxP - ' ^^y ■ (A. 6a) 

M = I^q - ‘ ^x^’^P (A. 6b) 

N = I^f - I^^(p - qr) - (I^ - Iy)pq (A. 6c) 

Kinematic Equations in Wind Axes : 

♦w " Pw *w ®w Pw '*’w ®W 

®w " % *w - I" '»w 
K " ^^’P ‘•'w Pw ®w 
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Without subscripts, the above equations also apply in body axes. 
Additional Kinematic Rel 


a = q - q^ sec 3 - p cos a tan 3 - r sin a tan 3 

(A. 8a) 

3 = r + p sin a - r cos a 
w ^ 

(A. 8b) 

p = p cos a cos 3 + (q - a)sin 3 + r sin a cos 3 
w 

(A. 8c) 

The velocity components relative to the earth fixed 
in terms of V are: 

reference system F^ 

Xr = V cos 0 cos ip + W 
t w w xL 

(A. 9a) 

= V cos sin 

(A. 9b) 

"e = -V ®w ^ ^zE 

(A.9c) 


For the body frame of reference, we obtain: 

= u cos e cos ij; + v(sin sin e cos \{j - cos (p sin i|j.) 


+ w(cos (p sin 9 cos \p + sin <p sin ijj) + W^^ (A. 10a) 

y^ = u cos 0 sin \{j + v(sin 4> sin e sin ip t cos (p cos 

+ w(cos (p sin 0 sin - sin 4 > cos ip) + (A. 10b) 

= -u sin 0 + V sin ij; cos 0 + w cos (p cos 0 + W^^ (A. 10c) 

Finally, the relationship between the velocity in the body frame of 
reference and that in the wind frame of reference is given by: 


u = V cos a cos 3 
V = V sin 3 
w = V sin a cos 3 


(A. 11a) 
(A. 11b) 
(A. 11c) 


Small Disturbance Theory with Variable Wind Field 

In most of the conventional analyses of aircraft motion, a linear- 
ized form of the equations, for small disturbances about a reference 
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condition of steady rectilinear flight over a flat earth, is employed. 
(Symmetric flight requires 7 to lie in the plane of symmetry). The 
linearized equations are developed by conventional methods; however, it 
will become apparent in the development that the reference conditions in 
the presence of a wind field are developed by conventional methods; 
however, it will become apparent in the development that the reference 
conditions in the presence of a wind field are difficult to define. The 
frame of reference for the small disturbance model is generally taken as 
the "stability" frame, with a special set of body axes coinciding with 
the wind axes in the reference condition, but departing from it and 
moving with the body during a disturbance. 

The steady state values of the variables are denoted by a subscript 
e, and changes from the steady state values are denoted by the prefix A, 
i .e. , 

V = + AV 

e 

(J) = (j)^ + A(j) (A . 1 2) 


etc. 

In this reference frame the state variables are normally taken as V , 

e 

a , 0 and il^ All other variables are zero in the reference state 
e we we 

and for these the prefix a is dropped. 

The small disturbance equations are now developed following Etkin 
[4]. The angle of climb is denoted by y, a more commonly used 
symbol. The angle is set equal to zero since initial heading has 
no special significance in the flat-earth approximation. This does not 
preclude the possibility of winds other than head-on wind, however, 
since the angle of the wind relative to the flight path is determined 
by the three components of the wind field. The thrust vector, T, is 
permitted to be at large angles to the direction of motion but is 
required to rotate rigidly with the vehicle when the vehicle is 
perturbed. Thus, in body axes; 
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fcos OLj) 


(T + AT) 


0 


^sin a^J 

and in wind axes: 

T^w = (T + AT) (cos 

% 

Ty^ = (T + AT) (-cos 
'*'zw " AT) (-cos 


cos a cos 3 + sin sin a cos 3) 
cos a sin 3 - sin sin a sin 3) 
sin a + sin a-j- cos a) 


(A. 13a) 
(A. 13b) 
(A. 13c) 


In the stability reference frame, = 0, hence, a = Aa. If* one makes 
the approximation sin A ~ A, cos A = 1.0 and neglects the squares and 
products of the a terms, Equation A. 2a becomes: 

(T + aT)cos - AaTg sin ay - D - AD - mg sin(Yg + Ay) 


= m(V + 1^1 ) + m(q W - r W ) 

xw^ zw w yw^ 

where the reference state is defined by: 

Tg cos ay - D - mg sin Yg = 0 


(A. 14) 


(A. 15) 


Under the assumption of uniform wind, = 0, the small disturbance 
approach is justified. On the other hand, for nonuniform wind fields: 


xw 


3W 1 

- + (V„ + AV) ^ 


3t 


3X 


(A. 16) 


w 


in the wind frame of reference, and 


'xB 


3W^ 3W 3W 

— + — u + — V + — — w 

3t 9X ^ 3y 3Z 


(A. 17) 


in the body frame of reference. Thus, a problem is encountered with the 
method of small disturbances for the case of a general wind field since 
a continual departure from the reference state with time occurs. 

If the wind is considered time dependent and the reference state is 
allowed to vary with time, then from Equation A. 2 the governing equations 
of the reference state become: 
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cos Ot - D - mg sin Y - nil^ = 0 


XW 


Cg + -"V = 0 


T sin + L + mg cos y + mlJ = 0 

t- I c G 


(A. 18a) 
(A. 18b) 
(A. 18c) 


These equations could be solved for y , C , and L., given a specified 

G 0 G 

wind field. The small disturbance equations for this time-dependent 
reference state then become: 

AT cos a-r - AaT^ sin a-j- - AD - mg Ay COS y = mV + m(W q 

I e I ^ ' e ^ zw^w 


- W , r ) 
yw Vi' 


-3T cos - aC + mgi , cos y^ = m[r, ,(V + W ) - p W ] 
e T ^^w 'e w' e xw' zvi 

AT sin aj + AaTg cos + aL + mg Ay sin y^ = - qw(V^ 


+ W )] 

xvi' 


{A. 19a) 
(A. 19b) 

(A. 19c) 


It is apparent, however, that the advantage of the small disturbance 
equations, which is that they are a linear time-invariant system of 
equations that can be solved by established mathematical transfer 
function techniques, is lost since the coefficients containing y^ are 
functions of time. 

A similar result is obtained with Equation A. 2, which with the 
small disturbance approximation, becomes: 


Ax - mg cos Ay cos Yg = mu + m(qW^ - rW^) 

Ay - mgtj) cos y^ = mv + m[r(Vg + W^) - pW^] 
Az - mgAy sin y^ = mw + m[pWy - q(V^ + W^)] 
where the reference state is such that: 


X =mgsiny -mW =0 
e ^ ex 


yg = mw^ = 0 


Zg + mg cos y^ - mW^ = 0 


(A. 20a) 
(A. 20b) 
(A. 20c) 

(A. 21a) 
(A. 21b) 
(A. 21c) 
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Linearized Equation of Motion for Uniform Wind 


For a uniform wind, = 0, and Equations A. 6 through A. 8 

A y z 


become: 

= ^xP - 

AM = I^q 
AN = I -r - 


\> = p + r tan y 


Y = q 


w 


ip = r sec Y 
w w ’e 


(A. 22a) 
(A. 22b) 
(A. 22c) 
(A. 23a) 
{A. 23b) 
(A. 23c) 


Without the subscript w, these equations apply in body coordinates. 
Also, we have the kinematic relationships: 


a = q - q 


w 


3 = r„ - r 


P„ = P - 


The aircraft velocity in earth coordinates becomes: 

Xg = Vg(cOS Yg - Ay sin Yg) + AV COS Yg + W^j. 

= Vg cos Yg + Wy£ 

Zr = -V sin Y« - aV sin y^ - Ay cos y^ + W^r 
E 'e e e e zE 

or in body coordinates: 

Xr = V cos y + V Ay sin Yo cos y« + w sin y^ + W r 
yn = Ve ♦ cos Yg + v + W^^ 

^E ■ ~^e ''^e ■ '^e Yg - u sin Yg + w cos Yg + 


(A. 24a) 
(A. 24b) 
(A. 24c) 

(A. 25a) 
(A. 25b) 
(A. 25c) 

(A. 26a) 
(A. 26b) 
(A. 26c) 


Recall that u^ = and = 0 in the reference state. The rela- 

e e e e 

tionships among the wind components in the earth frame of reference and 
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the wind frame of reference become: 


W 


xw 


= W c(cos Y - Ay sin y^) + W 


'yE 


cos Yg “ W^^(sin Yg 


W. 


yw 


zw 


+ Ay cos Yg) 

= "''' ^e - V •' ^zE'^w ^e 

= W^^{sin Yg + Ay cos Yg) + ^e " 

+ ^^^(cos Yg - Ay sin Yg) 


(A. 27a) 
(A. 27b) 


(A. 27c) 


The equations are valid for the body frame of reference without the 
subscript w. 

Conventionally (i.e,, without atmospheric motion, 0=0), Equations 
A. 19a, A. 19b, A. 21b, A. 23b, A. 24a, A. 25a, and A. 25b are taken to be the 
longitudinal equations since they contain only longitudinal variables 
(aV, Aa, q, Ar, x^, z^) and Equations A. 20b, A. 21a, A. 21c, A. 23a, A. 23c, 
and A. 26b are taken to be the lateral equations since they contain only 
lateral variables (v, p, r, <p, ip, y^) . The equations thus decouple and 
form two independent sets which can be solved separately. However, with 
a wind, even a uniform wind, the longitudinal equations do not separate 
because p^ and r^ appear in Equations A. 19c and A. 19a, respectively. On 
the other hand, the lateral equations separate in view of the fact that 
neither 

'*■ *^x^ " '^xE "^e " *'^zE ^e^ (A. 28a) 


nor 

pW^ = "^e ^zE "^e^ (A. 28b) 


contain any of the longitudinal variables. 

Finally, the special case of a horizontal wind oriented parallel to 
the direction of motion, i.e.. 


fw 


t 


xE 

0 

0 


(A. 29) 
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results in a form of the equations which permits separation of the longi- 
tudinal equations as well. The equation thus has the familiar form: 


AT cos Yy - AaTg sin a-j. - AD - mg Ay cos y^ = mV + niW^^q^ sin y^ 
AT sin + AaTg cos tty + aL + mg sin Yg = cos Yg) 
AY + mg(J) cos Yg =" niv + m[r(Vg cos Yg) - sin Y^] 


(A. 30a) 
(A. 30b) 
(A. 30c) 
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APPENDIX B 


NOMENCLATURE 

\ap Acceleration capability 

C Side force 

D Drag force 

DA Acceleration difference 

F Frame of reference 

GNS Ground speed 

GS+ Flight deterioration parameter (Table 4.2, 2a) 

GS- Flight deterioration parameter (Table 4.2, 2b) 

g Gravity 

H altitude of airplane CG 

H6 Height of the glide slope above the ground 

HP Height of the aircraft above the ground 

h^ Arbitrary reference height scale 

I Moment and/or product of inertia 

IAS Indicated airspeed 

L Lift force 

L Rolling moment 

L Monin-Obukhov stability length scale 

L Turbulence length scale 

M Pitching moment 

m Mass 
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N Yawing moment 

p Rate of roll 

q Rate of pitch 

r Rate of yaw 

T Thrust 

T Time period 

TAS True airspeed 

u x-component of aircraft velocity relative to the atmosphere 

u* Friction velocity 

V+ Flight deterioration parameter (Table 4.2, 4a) 

V- Flight deterioration parameter (Table 4.2, 4b) 

V Relative velocity vector (airspeed) 

V Selected approach speed (kts) 

Airspeed 

V_ Approach airspeed 

®0 

V, Stall airspeed 

®stall 

— y 

Inertial velocity vector 

V y-component of aircraft velocity relative to the atmosphere 

iJ Wind velocity vector 

W Mean wind speed 

WD Difference in wind speed at the runway and at the aircraft 

WXgnd Wind component at the ground 

w z-component of aircraft velocity relative to the atmosphere 

X X-component of aerodynamic force 

X Distance along x-axis 

V Y-component of aerodynamic force 
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y Distance along y-ax1s 

Z Z-component of aerodynamic force 

z Distance along z-axis 

z^ Surface roughness 

Greek Symbols 
a Angle of attack 

3 Angle of yaw 

y Pitch angle 

r Wind shear vertical gradient in horizontal wind (aW /9z) 

X 

aA Acceleration margin 

Ax Deviation from desired touchdown 

9 Euler angle (elevation) 

K von Karman constant 

A Wavelength 

a Wind shear parameter (V (aW /az)/g) 

a Turbulence intensity 

<|) Euler angle (bank) 

i{j Euler angle (azimuth) 

ij;(z^/L) Stability parameter 

o)ph Phugoid frequency 

Subscripts 

i Initial value 

e Reference state 

E Measured in the inertial coordinates 

0 Landing speed 
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s 

T 

w 

X 

y 

z 


stall speed 

Direction of thrust 

Measured in the wind coordinates 

Measured in x-di recti on 

Measured in y-direction 

Measured in z-di recti on 


Superscript 


Time derivative d( )/dt 


Prefix 

A Small perturbation 
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of the effects of wind shear on aircraft performance are presented. These are 
followed by a review of mathematical solutions to both the linear and nonlinear 
forms of the governing equations. Solutions with and without control laws are 
presented. 

The application of detailed anlaysis to develop warning and detection systems based 
on Doppler radar measuring wind speed along the flight path is given. A number of 
flight path deterioration parameters are defined and evaluated. Comparison of 
computer-predicted flight paths with those measured in a manned flight simulator is 
made . 

The report ends with a review of some proposed airborne and ground-based wind shear 
hazard warning and detection systems. The advantages and disadvantages of both 
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